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SUMMARY 


A theoretical and experimental study was conducted to develop a validated 
first principles analysis for predicting noise generated by helicopter main-rotor 
shed vortices interacting with the tail rotor. A generalized noise prediction 
procedure was formulated in which the incident vortex velocity field, rotor 
geometry, and rotor operating conditions are input variables. The analysis includes 
compressibility effects, chordwise and spanwise noncompactness, and treats oblique 
intersections with the blade planform. Output from the analysis provides acoustic 
spectra and pressure signatures as a function of observer position. 

Assessment of the theory involved conducting a model rotor experiment which 
isolated the blade-vortex interaction noise from other rotor noise mechanisms. 

The experiment was conducted with a 0.76 m articulated rotor operating in the 
United Technologies Research Center Acoustic Research Tunnel. An isolated tip 
vortex, generated by an upstream semispan airfoil, was convected into the model 
rotor to generate blade-vortex interaction noise. Acoustic spectra, pressure 
signatures, and directivity were measured for different combinations of rotor tip 
Mach number, blade pitch, and blade number. Since input to the acoustic predic- 
tion required a knowledge of the vortex properties, blade vortex intersection 
angle and intersection station were determined using smoke visualization. Vortex 
strength and vortex core radius were documented with hot film probe measurements. 

Ingestion of the vortex by the rotor was experimentally observed to generate 
harmonic noise and impulsive waveforms. Assessment of the theory using the measured 
data verified the linear dependence on blade number and vortex strength. The 
absence of noise for vortex trajectories normal to the blade was also predicted. 

For oblique intersections with a blade, changes in acoustic radiation, due to inter- 
section angle modifications, were calculated to within 2 decibels of the measured 
data. A simplified scaling law was also developed to predict these changes. Low 
sensitivity to blade pitch angle modifications for typical rotor operating condi- 
tions was confirmed experimentally. The strong dependence on local Mach number 
was predicted although the sensitivity to Mach number changes was overpredicted at 
low frequencies. 

Measured blade vortex interaction directivity patterns confirmed the theoretical 
dipole acoustic source model including the location of the minimum in the plane of 
the blade. Absolute levels for acoustic spectra and directivity patterns were cal- 
culated, without the use of empirical or adjustable constants, to within 5 decibels 
at low frequencies. At high frequencies, sound pressure levels were underpredicted. 
General features of the measured vortex interaction acoustic pressure signature 
were calculated but the exact waveform shape and amplitude were not predicted. 



INTRODUCTION 


Background 

The noise level at a given observer position relative to a helicopter 
and within a specified frequency band represents the combined effect of a 
number of separate acoustic source mechanisms. Some mechanisms occur only 
on the main rotor, others apply to both the main rotor and tail rotor, and 
in more complex situations, interactions occur linking the wakes of one rotor 
with noise generated by the other rotor. The simultaneous occurrence of the 
acoustic source mechanisms has complicated the design of quiet helicopters . 
and contributed to empiricism and extrapolation of measured data when pre- 
dicting full scale helicopter noise signatures. 

Due to the limitations of empirically based helicopter noise prediction 
procedures, an increasingly important objective of helicopter noise research 
has been the development of procedures based on first principles analyses. 
This effort has been reasonably successful for a number of noise mechanisms 
such as blade thickness, loading, turbulence ingestion, and boundary layer- 
trailing edge interaction. However, in the area of main-rotor shed vortices 
interacting with the tail rotor, first principle analyses did not exist. The 
lack of a firm understanding of the tail-rotor vortex interaction mechanism 
created a need for an experimentally validated prediction procedure. In an 
effort to achieve this objective, the present joint experimental and theoret- 
ical research program was conducted to assess this noise mechanism. 

The blade-vortex interaction mechanism addressed here is considered to 
be potentially important. As illustrated in Figure la, under forward 
flight or vertical descent operating conditions vortices shed by the main 
rotor pass through the tail rotor resulting in the generation of noise. In 
principle, the strength of this noise mechanism can be reduced by changes in 
the helicopter design. For example, altering the relative position of the 
main rotor and tail rotor is a possible solution. In addition, changing the 
main-rotor blade tip shape may reduce the strength of the tip vortex or dif- 
fuse the vortex velocity field in a manner which reduces interaction noise. 

The decision to pursue such aerodynamic design changes is based on the 
ability to predict quantitatively the expected benefits. The major objective 
of the present study was, therefore, to develop and validate a first prin- 
ciples prediction procedure which starts with the incident vortex velocity 
field as an input. 
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Previous Investigations 


Experimental Studies - Few experimental studies of helicopter tail-rotor 
vortex ingestion noise are currently available. Leverton, et al. (ref. 1) 
observed that the Lynx helicopter produced such noise during approach condi- 
tions when the tail rotor intersected tip vortices shed by the main rotor. 

The experimental results indicated that this noise mechanism can be subjec- 
tively annoying and can dominate over the main rotor and engine noise. 

In addition to full scale measurements, Pegg and Shidler (ref. 2) studied 
the parameters controlling main-rotor tail- rotor interaction noise using 
a small scale helicopter model. The model permitted varying tail rotor posi- 
tion relative to the main rotor in addition to the tail rotor speed and rota- 
tion direction. Experimental results showed that certain features of the 
tail rotor acoustic spectrum were attributable to interactions with tip vor- 
tices shed from the main rotor. 

A limited study of model scale vortex interaction noise was also con- 
ducted by Paterson and Amiet (ref. 3). Interaction of an isolated tip vortex 
with the model rotor was observed to produce significant harmonic noise 
which extended to high frequency. The impulsive acoustic waveform generated 
by the interaction was found to be sensitive to vortex strength and position. 

Theoretical Studies - Analytical models predicting the blade vortex in- 
teraction noise mechanism presently treat only the main rotor noise source. 

To help differentiate these analyses from the present formulation, a brief 
discussion of the main rotor interaction mechanism will be given. In this 
case the blade and vortex lie in different (or identical) parallel planes.. 
Furthermore, the blade leading edge is parallel to the vortex centerline in 
the most intense interactions (Fig. lb). This condition is modeled by a two 
dimensional analysis assuming end effects can be neglected. On the other 
hand, if the vortex centerline is skewed relative to the airfoil leading edge 
(but still parallel to rotor plane) , the vortex velocity field sweeps across 
the airfoil. Under this condition, a three dimensional analysis is required. 


The above described two dimensional blade vortex interaction was treated 
by Amiet (ref. 4) assuming linearized flow. The three dimensional interaction 
was modeled by Widnall (ref. 5), Widnall and Wolf (ref. 6) and Adamczyk (ref. 
7) . These analyses include the two dimensional problem as a special case 
but they do not treat the present tail rotor interaction geometry which in- 
volves a nonplanar intersection. In this case, the rotor disk is approxi- 
mately at 90° to the vortex plane at the intersection station. One objec- 
tive of the present study was, therefore, to develop the specific analyses 
needed to treat the tail rotor blade vortex interaction geometry. 
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Present Investigation 


Ob.j ective - The overall objective of the present study was to develop 
an experimentally validated blade-vortex interaction analysis. The analyt- 
ical phase of the program required the formulation of a generalized noise 
prediction procedure in which the incident vortex velocity field, rotor 
geometry, rotor operating conditions, and far field observer position were 
input variables. Empirical scaling factors were avoided in this first prin- 
ciples analysis. 

Assessment of the theory required conducting a tractable acoustic exper- 
iment which isolated this noise mechanism from other operative rotor noise 
mechanisms. Also, detailed measurements of the incident vortex velocity 
field were obtained to provide the experimental input to the analysis. Con- 
firmation of the theory for the experimental test conditions involved assess- 
ment of blade vortex interaction noise sensitivity to changes in the vortex 
field and the rotor operating conditions. 

Theoretical Approach - The present analysis transforms the incident vor- 
tex velocity field into a spectrum of sinusoidal gusts incident on an isolated 
rotor blade. Linearized theory is then used to determine the unsteady blade 
lift response function. The response is given by a generalized function 
which includes compressibility effects, chordwise and spanwise noncompact- 
ness, and treats skewed or normal gusts in the decomposition of the incident 
vortex field. Input to the analysis requires a knowledge of the vortex 
strength, vortex viscous core radius, blade-vortex intersection angle, blade 
vortex intersection station, blade tip Mach number, blade pitch angle, blade 
chord, blade number, and far field observer location. Output from the analy- 
sis provides acoustic spectra and pressure signatures as a function 
of observer position. 

Experimental Approach - Figures 2 and 3 show the experimental arrange- 
ment employed to assess the present analysis. Here a tip vortex, generated 
by an upstream semispan airfoil model, is convected into a model helicopter 
rotor operating in vertical ascent with a 9 m/sec inflow velocity. Although 
this test configuration lacks the appearance of a tail rotor, it represents 
a tractable geometry which permits direct assessment of the blade-vortex in- 
teraction mechanism. Selection of this geometry was based on an evaluation 
of the parameters controlling the noise mechanism and a detailed discussion 
is included in the description of the experiment. 

The experiment was conducted with a 0.76 m diameter articulated rotor. 
Acoustic spectra, pressure signatures and directivity were measured for dif- 
ferent combinations of rotor tip Mach number, blade pitch, and blade number. 
Blade vortex intersection angle and blade intersection station were documented 
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with smoke visualization. Vortex strength and core radius in addition to 
viscous core turbulence intensity and length scale were documented with hot 
film probe measurements. 

It was considered important to devise an experiment in which only one 
parameter varied at a time. This was necessary since the vortex and rotor 
aerodynamic parameters were coupled. For example, changing the rotor tip 
Mach number also changed the mean flow contraction ratio for the rotor orien- 
tation shown in Figure 2. Hence, the incident vortex strength, intersection 
angle, and intersection station also changed. Selection of the test matrix, 
therefore, required documentation of these coupling effects. 
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LIST OF SYMBOLS AND NOMENCLATURE 


a 

b 

B 

BPF 

c 

c o 

c 2 

D 


E* 


F 

f 

f(M) 


G,G 


angle between hot film sensor and total velocity vector 
semichord 

rotor blade number 
blade passing frequency 
chord 

sound speed 

upstream semispar airfoil chord length, 11.4 cm 

vortex core diameter 

combination of Fresnel integrals 

point force 

frequency, cycles/sec 

function of Mach number; see Eq . (13) 
normalized pressure jump; see Eq. (2) 

single sided spectrum for far field sound generated by single 
and multiple vortex intersections 


H vertical separation distance between tip of semispan airfoil 

and tip of rotor blade 

h semispan airfoil open jet penetration height 

I distance from rotor blade tip to vortex intersection station 

i v/=T 


J Y 
n n 

J 

ux 

K 


K 


1 


Bessel functions of the first and second kind 
Eulerian time scale in autocorrelation 

coefficient relating velocity and voltage for hot film probe 
measurements 

vortex strength 
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K 


2 


K 

y 


k 

X 



k * 

X 


k 

x°° 

k 

k 

JC 

M 

Mr 

M 

00 

p 

AP 

R 


R 1 

RPS 


r A»' r B 


S 


normalized vortex strength, = K^/cU o 

2 

specific value of k^; uS y/cr 
chordwise and spanwise gust wavenumbers 




2 2 
k 3 /6 


x 00 

magnitude of wavevector 

unit vector in z direction 

effective lift function; see Eq. (6) 

local Mach number 

rotor blade tip Mach number 

similarity Mach number for skewed gust 

pressure 

pressure jump across airfoil 
observer-vortex distance 

effective radius of exponential decay function applied to vortex 
potential field 

rotor radius, 0.38 m 

autocorrelation functions for far field sound generated by a 
single and multiple vortex intersection 

rotor revolutions per second 

distance from observer to airfoil-vortex intersection point. 
Also, local radial station in vortex 

radius of vortex viscous core 

radius of open jet flow before (r^) and after (r^) contraction 
induced by rotor 

Sears function; see Eq. (13) 
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s, ,s 

1 cr 


SPL 

T 

t 

U 


U N1’ U N2 


V 


VVV V 4 


V 


N 


V, 


v 


w 


w 


x,y,z 


x 


v y i 


two sided spectrum for far field sound generated by a single 
and multiple vortex intersection 

2 

sound pressure level in dB relative to 0.0002 dyne/ cm 
time increment in Fourier summation 
time 

fluid velocity relative to airfoil in present theory; c Q M . 

Also, axial velocity in vortex core or in semispan airfoil two 
dimensional wake 

velocity components normal to hot film probe sensing element in 
two different orientations 

total velocity vector in the plane of the hot film sensor 
sweep speed of gust along airfoil relative to airfoil 
vortex azimuthal velocity vectors 

vortex velocity component normal to blade planform 
sweep speed Of gust along airfoil relative to fluid 
vortex azimuthal velocity 
maximum vortex velocity 

sinusoidal gust velocity normal to plane of airfoil 
Fourier amplitude of w^ 

Cartesian coordinates; see Fig. 4 
airfoil coordinates 

angle between gust and airfoil. Also, constant used in Eq. (22) 
rotor blade pitch angle 
semispan airfoil angle of attack 
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1 


r 

y 


6 


0 

1 


®2 


0 

o 


e 

V 




y 


V 


£,n,c 


p 

o 

a 

4> 


(0 

$2 

Overbar 


experimentally determined blade vortex intersection angle in 
vertical plane 

angle between hot film support prongs and sensor 

Vl-M 2 

OO 

circulation 

experimentally determined blade vortex intersection angle in 
horizontal plane 

Dirac delta function 

parameter defined by Eq. (20a) 

parameter defined by Eq. (20b) 

observer angle; see Fig. 4 

vortex angle; see Fig. 4 

Mk /e 2 
x 

M k /I3 2 

CO X 

Hk /gk 
x y 

Cartesian coordinates; see Eq. (37) 
freestream density 



velocity potential 
observer angle; see Fig. 4 
vortex angle; see Fig. 4 
radian frequency 

radian frequency increment in Fourier summation 
quantity normalized by the semichord b 
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ANALYTICAL FORMULATION OF THE BLADE VORTEX INTERACTION MECHANISM 


Definition of the Problem and Description of Approach 

The present analysis was developed for a blade vortex interaction 
geometry in which the rotor plane of rotation and the vortex plane are 
skewed. Consequently, the blade cuts through the vortex centerline during 
its motion as illustrated by the combined airfoil and vortex geometry shown 
in Figure 4. The orientation of the vortex relative to the airfoil is arbi- 
trary. Small interaction angles between the airfoil plane and the vortex 
centerline are, however, not permitted in the analysis. This restriction 
occurs because the intersection point of the airfoil and the vortex is 
assumed to be localized on the airfoil allowing end effects to be neglected. 
Thus, the blade vortex orientation treated in references 1-4 is excluded from 
the present analysis. 

The present formulation is based on linearized theory with the rotor 
blade modeled as a flat plate of infinite span. Because of the linearized 
flow assumption, steady loading is assumed to be small and is not coupled 
with the unsteady blade vortex interaction noise mechanism. Standard airfoil 
theory dictates that only the fluid velocity normal to the airfoil plane 
affects the airfoil loading. Also, because of the linearized flow assump- 
tion, the vortex velocity field and its trajectory are not altered by velocity 
perturbations induced by the airfoil. 

Analytical expressions used for the airfoil response functions are for 
the case of a flat plate airfoil of infinite span. A sinusoidal gust upwash 
is assumed incident on the airfoil. Thus, the airfoil response functions are 
similar to the classical Sears function (ref. 8). The only difference between 
the airfoil response functions used here and the Sears function is that the 
present solutions extend the Sears solution to the case of compressible flow. 
Also, using the similarity rule of Graham (ref. 9), the gust can be skewed 
relative to the airfoil. If the Mach number is set equal to zero and the 
gust is assumed parallel to the airfoil, the present solution reduces exactly 
to that of Sears. 

A detailed derivation of the present gust response functions is given 
in references 10-14. Separate solutions are given for the low and high fre- 
quency regime. The low frequency solution is a series type of solution and 
the high frequency solution is based on a convergent iteration procedure; the 
solution produced is asymptotically correct as frequency goes to infinity. 
Comparison with exact numerical solutions of Graham (refs. 9, 15) shows that 
the present solutions give an error in the lift of at most a few percent. 

The maximum error occurs near the changeover point from the low to the high 
frequency solutions. 
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The model for the vortex velocity field is the same one used by several 
other authors (refs. 16, 17). It reduces to a potential vortex at large dis- 
tance and approaches a solid body rotation in the core. This avoids the in- 
finite velocity existing at r = 0 in a potential vortex. 


General Expression For Far-Field Sound 

The response of the airfoil to the vortex velocity field will be treated 
by Fourier decomposing the velocity field and calculating the airfoil response 
to each resulting gust component. The inverse Fourier transform then gives 
the time behavior of the gust-vortex interaction. 

The airfoil is assumed to be a flat plate of infinite span lying in the 
x,y plane, with the leading edge at x = -b (see Fig. 1). If the airfoil 
lift were being calculated, the assumption of infinite span would create prob- 
lems; since the vortex field has only an r“^ decay rate, an infinite lift 
would be obtained in general. This singularity would become apparent on 
attempting to perform the Fourier inversion of the expression for lift pro- 
duced by a gust component; this has a k~^ wavenumber behavior for small k. 
However, since the present analysis calculates the sound and the sound is 
dependent on the time derivative of the loading, an additional factor of k 
is introduced which eliminates the singularity at k = 0. There remains, 
however, a difficulty in the calculated value of the average pressure level. 
This is considered in the section titled Experimental Assessment of Blade 
Vortex Interaction Noise. 

Linearized flow is assumed. This means that the vortex remains fixed 
relative to the mean flow at infinity. Thus, on Fourier decomposition the 
z component of the vortex velocity field in the x,y plane can be written in 
terms of the gusts 

w g = w(k x ,ky) e "'( k x( x_ut ) + k y y j 

This gust will lead to a pressure jump across the airfoil which can 
be written 


AP( x ,,y | ,k x ,k y )= 277 y> 0 Uw(k x , 


k y )g(X|,k xl k y )e- ,{k yyi^xUi) 


( 2 ) 


The far-field pressure of a tranverse point dipole or force F k e at 
x l> yi is 


P = — ■■ Fhiz - e '[<«t+M(MX-0-)J e ->MfM*Hxx,+/ 3 !! yy | )/ 0 -] 

4ttc 0 ct z 1 


(3) 
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This is obtained from the usual expression for a dipole (ref. 11) by 
neglecting terms OCr^/r) where r^ is a source (airfoil) dimension and r is 
the source-observer distance. 


The far-field pressure due to the gust component w g is found by replac- 
ing F in Equation 3 with AP from Equation 2 noting that 

CJ=k x U ( 4 ) 

and integrating over the airfoil surface. This gives 


P(k x ,k y )= 


i77 7>o u2bk x z 


c o CT 


8 (k y -^.£ 2 y/cr) w(k x ,k y )/(k x ,k y ,M) e 


i k x Ut+fj.(Mx-<r)](5) 


where 

./(k x , k y , M) = -t f gU, ,k x ,k y M) dx, (6) 

represents an effective lift per unit span. The exponential phase factor 
in £. accounts for retarded time effects for sound propagating from source 
to observer; if this factor were set equal to 1, Jl would represent the actual 
lift per unit span. 


To find the far-field pressure-time response, the integral over all k. 


and ky gust components is taken giving 


P(t) = - 7 T 


5 £^_/®[- i k x S(k x ,k y )]j(k xl K y ,M). , [ k * U,+ ^ ,,,l -' , l dk] 


(7) 


where 

K y = ^/3 2 y/o- 
fi*k x M//3* 

/3 Z = i - M 2 

cr 2 = x z +/3 2 (y 2 + 2 Z ) 


2 

The factor exp [iy (Mx-a) ] gives only a time shift At = (a-Mx)/3 c to account 
for the time for the sound to propagate from the airfoil leading edge to the 
far field and can be ignored. Thus, 


P(t) g — [-ik x w(k x ,K y )]^(k x ,K 


-'to 


M )e ik ^ 


dk. 


( 8 ) 
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The Airfoil Lift Response 


Graham (ref. 9) has shown that for a flat plate infinite span airfoil in 
subsonic flow the response function to a skewed gust can be obtained by similarity 
considerations from either (1) the gust response of the airfoil to a parallel gust 
in compressible flow or (2) the gust response of the airfoil to a skewed gust in 
incompressible flow. The first or second solution should be used depending on 
whether the intersection point of the gust wavefront with the airfoil leading edge 
moves supersonically or subsonically relative to the fluid, respectively. 

For the present case it can readily be shown that the intersection point 
always moves supersonically for the gusts of interest. Basically this is because 
a subsonic gust leads to an exponentially decaying pressure field away from the 
airfoil. The angle a of the gust wavefront relative to the airfoil leading edge 
is 


tan a = Ky /k x = My /cr 

The intersection point sweeps along the airfoil with velocity 


(9) 


V =u cot a = c 0 cr/y 

Relative to the fluid the velocity is 


( 10 ) 


v r = J '/u 2 +V 2 = c o J x 2 + y 2 +/3 2 Z 2 /y - C 0 


( 11 ) 


Thus, the solution requires only the airfoil response to a parallel gust in subsonic 
flow, along with the similarity rules of Graham. 


The similarity rule of Graham (ref. 9) relates the surface pressure gust 
response functions g(£, k x , ky, M) for a skewed gust in compressible flow to the 
function g(£, k x , 0, M) for a parallel gust in a flow of Mach number M. The rela- 
tion is 


g(£,k x ,k y ,M) = 


“gie.k 


X GO > ^ CO ^ e 


•£k v 2 /k. 


( 12 ) 


where 

k * go s k x /3 cq 2 //3 2 



I- M 


2 

CD 


M C0 = M V j - V “ 2 
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* = Mk x //3 k y 

^CO ~ ^ CO 00 ^00 " ^ GO ^ x /ft 


The response function for a parallel gust is divided into two ranges in which 
different approximate solutions are used. 

Case 1 : p,*, - 0.4 For low frequency an approximate gust response function is 

(ref. 10) 


where 


g(x,k x ,o,M) = 



s(k *)e 


ik, 


[M 2 x + f(M)] 


k x* = k x //3 2 

f (M) — ( I - fi) M + {3 (l + j3) - 2 


(13) 


and the bar over x and k x indicates normalization by the semichord b. The function 
S is the Sears function (ref. 8) used in incompressible flow airfoil theory. The 
above solution ignores terms 0(M k x ) but retains lower order terms such as 0(k x ), 
0(M 2 ), 0(k* M^) . Substituting Equation 13 for g(£, k xoo , 0, Moo) in Equation 12 gives 


g(x,k x ,K y ,M)= s(k*)e ifiMX+ ' k * f(M “ 1 


(14) 


Together with Equation 6 this gives for J. 


^(k x , K y ,M) = -±-5(k *)e * 


(15) 


Case 2: \i w > 0.4 For high frequency the airfoil response can be calculated 

using a method of Landahl (ref. 19) which treats the airfoil as a semiinfinite flat 
plate with a leading edge but no trailing edge. A correction can be added which 
treats the airfoil as a semiinfinite flat plate with a trailing edge but no leading 
edge. An iteration can be carried out between the leading and trailing edges and 
was shown by Landahl to converge for all frequency. The first two terms of this 
solution are (refs. 7, 11, 12) 


g,(x,k ,o t M) = 


■nr y^r+My^ri ) 


i[/2(!-M)(l + x) + 7r/4-k x ] 


(16a) 
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g 2 (x,k o,m)« 


tt \/2tt(| + M)k 


— [(I +i) E*(2^i(l-x))-|] e _i [^* l_M )( | +*)+TT/4-k x ] 


(16b) 


where 


E* = ! 


%-< 1L. 

JTn ■£, vT 


/ 


(17) 


E* is a combination of Fresnel integrals. Substitution in Equation 12 gives 
for a skewed gust 

9 |(x,k x K yi M) = 


7 t rrk x (l + M co )(l + x) 


g 2 (x,k x ,K yl M) ~ 


7T *y 27rk (1 + Mqq) 


ik x * [mJi+x) - M 2 x- |]- i tt/ 4 
s -ik x * + -tj-iir/4 


(18a) 


[(l + i) E*(2k x *M c0 (l-x))-l] 


(18b) 


From Equation 6 Jl is found to be 

^,(k xt K y ,M} = 

^ 2 (k x ,K y ,M) « 


E*(2©,)e i02 


77" v^k x (H-M 00 )©, 


(19a) 


'©2 


7 T ©^TTkj^d + M^j) 

r * - * yF e~ i2 ®i 

+ (i-i)|E*l4k x *M 00 )- 


{i(!-e"' 201 ) 


where 


1 +(M/M 00 )(x/cr) 




(19b) 


r\ - / M C0 X \ 


(20a) 


©2= k x + (M- x/cr) - 7T/4 


(20b) 
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The above equations for g and I are for an airfoil situated between -b<x<b 
whereas the computer program calculates values for £. for an airfoil situated between 
0<x<c. 


The Velocity Upwash 

The velocity field produced by an ideal vortex in an incompressible fluid is 


_ r 

V 8 2irr 


( 21 ) 


where T is the circulation and T / 2 tt is the vortex strength. A more realistic model 
which includes the effects of vorticity and avoids the singularity at r = 0 (by 
introducing a decaying radial function to simulate the viscous core) is 


V 
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_ e - Q<r/r o> 2 ] 


( 22 ) 


where a = 1.25643 is determined by requiring dvg/dr = 0 at r = r Q . Also, v^ = v Q 
for r = r Q and the circulation T m for large r is r w = 2tt(1 + 0.5/a)r o v Q . This vor- 
tex model has also been used by other authors (refs. 16, 17). 

For a flat plate airfoil, only the gust velocity normal to the airfoil surface 
produces a lift response. If the airfoil were to intersect the vortex normal to 
the vortex line, no lift would be produced. Thus, it is necessary for the vortex 
to be tilted out of the normal. The vortex will be assumed to lie in the x,z 
plane. This will later be extended to permit the vortex to be rotated by the polar 
angle <f> v about the z axis as well, with the airfoil span along the y axis. 

To find the z velocity component for a vortex tilted by 0 V from the z axis and 
4> v = 0, the radius r in Equation 22 must be interpreted as the normal distance R to 
the vortex; i. e. , 


R 2 =(x cos 0 v - z sin 0 v ) 2 + y 2 


The velocity is then 


(23) 



V R > 


(24) 
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For a potential vortex Equation 21 together with Equation 24 give 


yr 

W P 0 ' ' 2ttR 2 


sin © v 


( 25 ) 


This can be readily checked by noting that the velocity potential $ for a potential 
vortex is 


<*> = 


2 IT 




(26) 


where $ represents the azimuthal angle around the z’ axis in a coordinate system 
with z' along the vortex. The relation between the primed and unprimed coordinates 


is 


Then, since v = VC>, 


x' = x cos 0 V ~ z sin 0 V 

y / = y 

z' = x sin 0 V + z cos 
ton <£ 0 = y W 


(27) 


<3 z 

agreeing with Equation 25 


. _E_ r iio Hi + Mo i*Li,_r_,ine 

27r L a 2 ' az ax' dz J zir 


v x' 2 +y' 2 


(28) 


In order to represent the upwash field as sinusoidal gusts, the x,y Fourier 
transform of Equation 24 will be performed. Thus, 

sin t 7 v JJ _y V 0 (R)e _i(k * x + k y y ) dxdy 


S (L k ) = — — 

VR *’ y' (2 tt) 2 1 


CO 


(29) 


Transforming from an x,y to an x',y' integration (noting that the integration is 
performed in the z - 0 plane) then to an R, <p integration gives 


w(K x ,ky)» t ° t1 | 2 v f C0 Rv e lFt)dF)f 2lr sin<t>e ,(k * cos8 v + k y s ' n< *’ ,R d <^> 

*27 rj 0 q 

The <p integral can be performed for arbitrary v 0 (R) . Denoting this integral 


(30) 


by I, 


1= Z\ sin sin[R(k x cos <£/ cos 0 V + k y sin <£)] d<£ 
o 

= 4i J cos(£ k x R/cos 0 V ) sin (R k y J I - £ 2 ) d £ 


(31) 
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This integral is shown on page 399 of reference 20 to be 


1 = 


27 rik y — __ J,(Ryk v 2 + k x 2 /cos 2 6 W ) 


Vky 2 + k x 2 /cos 2 6 


( 32 ) 


Substitution of this result into Equation 30 and using Equation 22 for Vq(R) 


gives 


w (k x , ky) = 


i ky $1 n Q y | 

2 .y^ + k/ccs^ ll+ ^ )roVo 


/“ [«-e- a(R/r o ) ] J,(Ryky Z + k x 2 /cos 2 ^ v )dR 


(33) 


Reference 20 on pages 717 and 940 gives the results 


/ CO ? i 

J, (/3 K)e~ a * dR«ir (1,^) 


(34) 


Equation 33 then becomes 

w(k x ,k y ) = 


2 

ro,— )=l-e-£ Z / 4a 

4a 


(35) 


r(i,co) = i 

ikyd+ 2 ^)t° n 0 V 


Zw (ky 2 + k x 2 /cos 2 0 V ) 


r 0 v 0 e'(ky 2 + k x 2 /cos 2 5v)r 0 2 /4a 


(36) 


This can now be introduced into Equation 8 which, together with the expressions for X 
given previously, determines P(t) . Note that $ varies linearly with changes in vor- 
tex strength (T/2 tt r 0 v 0 ). 


The Velocity Upwash for an Arbitrary Vortex Orientation 

The velocity upwash given by Equation 36 is for a line vortex in the x,z plane 
while the airfoil is along the y axis. It is straightforward to generalize this 
to the case of an arbitrary orientation of the vortex relative to the airfoil. 

Let £, q, £ represent a coordinate system rotated about the z axis by 4> v . 

Then coordinates x, y, z are related to £, X, by 

x = £ cos sin d) v 

(37) 

y = £ sin c£ v + 17 cos <£ v 
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(38a) 


The relation between the wavevectors k x , ky and k^, k^ for the two coordinate 
systems are easily related using the basic definition 

0 k £ , k ) = (2V) 2 // w (C ,i») e- i(k ^ + S T)) d$ dT, 

' ^ -00 

w (k x ,k y ) = Jf f T 7 )e" ,(kxX+kyy) dxdy 

-co 

= (— ) 2 f /* w(C,i7)e' i ^ (k>COS * v + k > Sin ^ vl + i,,(kxSin ^ v_k > COSev) dCd'rj 
^ -00 

where the double bar over w denotes a Fourier transform in C and n and the 
double wavy bar indicates a Fourier transform in x and y. 

Comparing these two equations shows that 


w (k^ .k ^) 1 w(k x ,k y ) 


(39) 


k^ = k x cos<£ v + k y sin v 
k 7 ? =-k x sin^ v + k y cos c£ v 


( 40 ) 


If k x and k in Equation 36 are replaced by k^ and k n , then Equation 40 will give 
w(k x ,ky) for an arbitrary vortex orientation. 


Application of the Theory 

Acoustic Signature Calculation - Prediction of the time dependent far field 
acoustic pressure signature generated by the blade vortex interaction mechanism 
was based on Equation 8. The integral in this equation was evaluated using a 
Fast Fourier Transform. The effective lift function, in this expression was 
determined from either Equations 15 or 19 depending on the wavenumber, k x . Wave- 
number spectrum, was based on Equation 36 with k x and ky replaced by k^ and k^ 
(eq. 40) in the case of an arbitrary orientation between the vortex and the 
airfoil. 

Acoustic Spectrum Calculation - The harmonic level of the tones generated 
by multiple chopping of the vortex was predicted using Equation 8. This was 
accomplished by recognizing that by dropping the integral over (k„ = w/u) 
there remains a function P^(cj) which is related to the sound spectrum of a 
single vortex intersection. Since P-^(w) represents a single blade vortex 
interaction, this function is continuous in the frequency domain rather than 
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representing a discrete function with non-zero values only at the blade passage 

r\j 

frequency and its harmonics. The continuity of P]_(w) allows extending the 
analysis to the case of multiple chopping of the vortex. This is accomplished 
by assuming that the vortex intersections are independent events. Then, if 
Rl(x) represents the autocorrelation of the far field sound of a single vortex 
interaction and R^t) the autocorrelation of an infinite series of intersections, 
R^ can be defined by an infinite series of the functions R^ separated by the 
time of blade passage; i.e.. 


CO 

R 00 ( = Z R.(r- 2 tt n/BXi) (41) 

n =~oo 1 

where B is the blade number and ft the rotational frequency in radians /sec. Taking 
the Fourier transform gives sound spectrum, S^, for multiple intersection 

CO 

S 00 (w)=BfiZ S,(nBXl)8(u>- nBX2) (42) 

n= -go 

Here S^(to) is the spectrum for a single vortex interaction and can be related 
to the known function P^Cto) by expressing the autocorrelation as 


Ri (t) = lim 
T —00 


T 

J P, (t) p, (t-T)d t 

-T 


Taking the Fourier transform gives 


(43) 


Sj(w) = B XI P, (w) P*(w) 


(44) 


Then, combining Equations (42) and (44) , transforming to a single sided spectrum 
function, G, and expressing the results in terms of unit frequency (rather than 
unit radian frequency) gives 


6(2 (f) = 2B 2 X1 2 I pj (nBXl) P*(nBXl)8 (45) 

n =-go ' ' 

This represents the final equation used to calculate the acoustic spectrum 
generated by multiple blade vortex intersections. Note that the delta function 
is 6(f) in Eq. (45) rather than 6 (w) in Eq. (42) and recall that 6(m) = 

6(2TTf) = 6(f) / 2 tt. 
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DESCRIPTION OF THE EXPERIMENT 


Acoustic Research Tunnel 

Operating Characteristics - The experimental study was conducted in the UTRC 
Acoustic Research Tunnel. A detailed description of the facility is given in ref- 
erence 24. The tunnel, shown schematically in Figure 5, is of an open-circuit, 
open- jet design. The inlet is provided with a high length-to-diameter ratio honey- 
comb section and a series of turbulence suppression screens. These features, in 
conjunction with a large tunnel contraction, provide a spatially uniform, temporally 
steady flow with a controlled test section turbulence level of approximately 0.2 
percent. Turbulence generators and grids can be inserted upstream of the nozzle to 
generate wake profiles and a range of turbulence levels in the test section. 

The open jet test section is surrounded by a sealed anechoic chamber 4.9 m high, 
5.5 m long (jet centerline direction), and 6.7 m wide. The chamber walls are lined 
with 0.5 m high fiberglass wedges which provide an anechoic acoustic environment 
above 175 Hz. Downstream of the test section the airflow enters a diffuser by way 
of a circular collector that has anechoic treatment on its flow impingement surface. 
The diffuser operates unstalled and thus is not a major source of background noise. 

To avoid tunnel fan noise from propagating upstream into the anechoic chamber, a 
Z-shaped muffling section with two right-angle bends and parallel treated baffles 
is located between the diffuser and the fan. The 1100 kW centrifugal fan exhausts 
to the atmosphere through an exhaust tower. 

Tunnel speed is determined from total pressure measurements upstream of the con- 
tracting inlet and static pressure measurements within the sealed anechoic chamber. 
Since losses are confined to the boundary layer, total pressure upstream and down- 
stream of the contraction are predicted and have been verified, to be equal. The 
test section velocity has been shown to be temporally steady. 

Open Jet Geometry and Test Section Arrangement - A circular inlet nozzle with a 
1.28 m diameter was used in the test program. Initial facility tests described in 
reference 24 identified an acoustic coupling between the inlet nozzle and the down- 
stream collector lip resulting in edge tones at high tunnel speeds. Triangular tabs 
were used previously to suppress the acoustic feedback mechanism. However, several 
facility changes performed in 1979 eliminated the acoustic coupling problem per- 
mitting the nozzle to be operated without tabs. 

Figure 7 shows the open jet test section arrangement. A NACA 0012 model airfoil 
was located at the exit plane of the open jet nozzle. The 11.4 cm chord airfoil 
generated a tip vortex which was convected into the downstream rotor. The model 
airfoil was mounted in a clamping device which permitted changing the angle of 
attack, <* 2 , and the vertical flow field penetration height, h. Variations in h 
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controlled the vortex radial intersection station on the rotor blade, (I), and the 
intersection angle, B-^. The plan view of Figure 7 shows a second vortex inter- 
section angle, y, which was also monitored. 

It was recognized during initial planning stages that changes in h and 

altered the tip vortex strength. The complicated coupling effects between h and 

a 2 in addition to the effect of changing rotor blade tip Mach number (M T ) and blade 
pitch angle, a^, are discussed in the subsection titled Test Program. 

The airfoil and rotor shown in Figure 7 were stationed on the open jet center- 
line with the rotor hub located 0.58 m downstream of the airfoil trailing edge. 

This station was selected based on the previous study of Paterson and Amiet 
(ref. 3) which employed the same rotor model. Preliminary measurements conducted 
during their investigation showed that the rotor ingested the open jet shear layer 
when stationed at approximately 0.9 m downstream of the open jet exit plane. 
Significant turbulence ingestion noise was generated under this condition. This 

problem was eliminated by moving the rotor upstream to 0.58 m (ref. 3). The ab- 

sence of turbulence ingestion - noise at this station was verified experimentally for 
the present tip Mach number operating range, which was only 10 percent higher 
than reference 3. 


Model Rotor 


The model helicopter rotor tested in the present study employed a 0.76 m diam- 
eter rotor. The articulated hub was provided with NACA 0012 untwisted blades of 
5.1 cm chord. Blade pitch angle was set with a blade profile template and an in- 
clinometer. Blade-to-blade variations in pitch angle were limited to - 10 minutes. 
This was sufficient to avoid extraneous tones at multiples of shaft frequency. 


Rotor tip speed (or revolutions per second, RPS) was set using an adjustable 
speed control with a 1 Hz resolution. A once per revolution short duration pulse 
generated by an optical sensor was used to monitor the speed. This time marker 
signal was also recorded simultaneously, on analog tape, with the acoustic signal 
from each microphone. The marker signal was later used as a trigger to initiate 
blade vortex interaction acoustic signature ensemble averaging. 


The small rotor blade chord tested here resulted in Reynolds numbers which 
varied, with rotor RPM, from 4 x 10 5 to 6.5 x 10 5 at the tip and 9.5 x 10^ to 
1*5 x 10 5 at the blade root. Over this Reynolds number range boundary layers 
on the blade remained laminar for the blade pitch angles investigated. Laminar 
flow over the blade trailing edge generated vortex shedding noise which was 
observed as a broad peak in preliminary acoustic spectra. The apparent broad- 
band nature of this noise mechanism was explained earlier by Paterson and Amiet 
(ref. 3), to be caused by the Doppler shift of discrete tones generated by blades 
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approaching and retreating from a fixed microphone. In addition, the spanwise 
varying velocity on the rotor blade generated a range of discrete tones which 
further broadened the spectrum. 

From a practical standpoint the vortex shedding mechanism has not been iden- 
tified as an important noise mechanism on full scale helicopters. Furthermore, 
preliminary acoustic measurements obtained during the present model study indicated 
that the broad peak masked the tones generated by the blade-vortex interaction 
phenomenon. It was, therefore, necessary to trip the boundary layers on the model 
rotor blades and force transition to turbulence to eliminate the laminar boundary 
layer vortex shedding noise mechanism. 

The boundary layer trip procedure employed here was identical to the pre- 
viously reported approach of Paterson and Amiet. The method represents a variant 
of the technique first reported by Hama in reference 25. In the present study a 
strip of adhesive aluminum tape was cut with pinking shears and installed on the 
pressure and suction surfaces of each blade at the 25 percent chord station 
(Fig. 3). Vortex shedding noise was eliminated as verified in the subsection 
titled Test Program. The present boundary layer trip approach has the advantage 
of being repeatable. It also avoids the question of uniformity on each blade and 
similarity between blades when using the distributed roughness approach. 


Instrumentation 

Acoustic Measurements - Far-field rotor noise was measured with 0.635 cm 
diameter microphones at grazing incidence. Frequency response on each microphone 
was flat over the 175 Hz to 10 kHz range of the present study. Microphone signals 
were amplified and then recorded on magnetic tape. Dynamic response for the FM 
tape recorder system was flat over the acoustic frequency range investigated. Cor- 
rections for atmospheric attenuation were not applied due to the small acoustic 
source-to-microphone propagation distance and the low acoustic source frequencies. 

Source directivity information was obtained by locating microphones on a 3 m 
or 2.5 m arc (Fig. 8a). The microphones formed a horizontal plane which passed 
through the rotor rig centerline. Measurement angles, 0 O , given in Figure 8a are 
referenced to the upstream rotor axis. Angles ranging from 0 Q = 60° to 140° pro- 
vided data from 30° forward of the rotor plane of rotation to 50° aft. 

A microphone was also located in a vertical plane to provide additional direc- 
tivity information. The unit was placed at 0 O = 130° and <f> 0 = 135 deg where <j> 0 
defines an angle measured relative to the horizontal plane through the rotor rig 
centerline (Fig. 8b) . 

It should be noted that the horizontal reference plane selected for 4> 0 
here is different from the reference plane used in the blade vortex interaction 
theory. In the latter case, the reference plane passes through the blade vortex 
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intersection station which is also indicated in Figure 8b. Since the vortex 
intersection occurred at approximately 0.32 m below the rotor centerline the 
difference between the two values of <f> Q is less than 6° for microphones on a 
3 m arc. Thus, the definition of the microphone orientation angle, <f> Q , can be 
considered to be similar for the present theory and experiment. 

It is important to recognize that the microphone measurement angle, 0 Q , in 
the experimental study coincides with the observer angle, 0 Q , in the theoretical 
formulation (Fig. 4) only if certain conditions are satisfied for blade pitch angle, 
ot^, and the azimuthal microphone angle, <j> . The conditions are aj = 0° and 
4>o = 0° or 180°. The latter requirement for <f> 0 was satisfied for the microphones 
located in the horizontal plane as indicated by the arguments in the above para- 
graph. But, varied from 5° to 14°. Values of 0 O associated with the experi- 
ment and theory, consequently, differed by a value of aj_. Input values of 0 O 
for the present theoretical prediction must, therefore, be increased by to 
account for this difference. In the general case where the above condition on 
<p Q is not satisfied, microphone measurement angles must be defined directly in 
terms of the blade referenced observer angles in Figure 4. Otherwise, an analy- 
tical transformation is needed to convert measurement angles to blade referenced 
angles . 

Spectrum analysis was conducted with an 800 line Fast Fourier Transform spec- 
trum analyzer (400 line in dual channel mode). The effective bandwidth was 12.5 
Hz for the 0-10 kHz analysis range. Bandwidths for other frequency ranges given 
in this report are determined by direct proportion. A high resolution capability 
in the frequency domain of the spectrum analyzer was also used to evaluate narrow- 
band random features of rotor tone noise at high blade passing harmonics. The 
analyzer bandwidth for this data is given in the corresponding figure. 

Vortex Velocity Field Measurements - Mean velocity measurements of the vortex 
velcoity field were obtained using a single0.025 mm diameter (with a 0.5 mm sens- 
ing length) linearized hot film probe with the sensor slated at 45° to the open 
‘jet centerline. The single element sensor was used to decompose the total vortex 
velocity vector into the azimuthal (Vq) and axial (U) components shown in 
Figure 9. The probe was oriented with the 0.5 mm sensing length parallel to the 
local Vq vector. Details of the voretx measurement procedure are given in 
Appendix A. It should be noted that a crossed probe could not be used for 
surveying the vortex since the large spacing between the two sensor elements 
(1 mm ) would average Vq information inside the small radius (0.4 cm) 
viscous core. 

Dependence of the vortex velocity field on radial distance from the vortex 
centerline was obtained by traversing the hot film probe through the vortex on a 
line perpendicular to the vortex center. These traverses were conducted using a 
two-directional linear traverse system (- 0.0013 cm position accuracy) coupled with 
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two variable speed stepping motors and a two digital encorders (0-0013 cm resolution). 
The traverse system, hot film probe and upsteram semispan airfoil are shown in 
Figure 10. 

Before conducting the above-described radial traverse through the vortex it was 
necessary to locate the vortex center. Although the approximate location of the 
tip vortex from the upstream airfoil was known, its precise position in space could 
not be determined in advance. Alignment of the hot film probe with the upstream 
airfoil tip was not an acceptable solution to this problem. This is because of the 
small size of the viscous core (approximately 0.4 cm radius) and the airfoil down- 
wash velocity. The latter phenomenon translated the vortex center as the flow was 
convected from the blade trailing edge to the rotor plane. A procedure was developed, 
therefore, to locate the vortex center at the rotor plane (Appendix A) . 

Turbulence intensity and length scale in the vortex viscous core were estimated 
from hot film crossed probe measurements. Limitations of the accuracy of this data 
due to sensor spacing must be recognized when interpreting this data. Linearizers 
and a sum and difference circuit were used in conjunction with a true root-mean- 
square meter to estimate the azimuthal and axial turbulence intensity in the viscous 
core. Similar to the mean velocity measurements, the crossed probe was traversed on 
a straight line through the vortex center to obtain the radial variations of the 
turbulence field. Turbulence integral length scale on the vortex centerline was 
determined using a 400 line real time correlator and probability analyzer. 

F low Visualization - The objective of the flow visualization phase of the 
study was to document the details of the blade vortex intersection. The vortex 
pathline at the rotor plane was determined by injecting smoke into the tip vortex 
at its inception point on the upstream semispan airfoil. Smoke was released from 
a small elongated slit on the rounded tip of the airfoil. Smoke ejection rate was 
adjusted to avoid generating a high speed jet which would alter the vortex formation. 
A closeup of the smoke ejection is shown in the photographs of Figure 11. 

The smoke pathline was documented with a camera mounted in the rotor plane 
of rotation. A short duration light pulse was used to illuminate the smoke. For 
each test condition photographs were obtained in the vertical plane to document 
the vortex intersection angle, 3^ , and the intersection station, I, shown in 
Figure 7b. In addition, photographs in a horizontal plane recorded the intersec- 
tion angle, y, shown in Figure 7a. This latter angle required documentation be- 
cause the strength of the semi-span airfoil downwash at the rotor plane was unknown. 

A general conclusion, however, was that y = 0° for all rotor operating conditions. 

All values of 3^ and y are based on averaging three photographs of each test condi- 
tion. Resolution was estimated at - 0.5°. 
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A sample photograph obtained in the vertical plane is shown in Figure 12. 

The upstream airfoil appears on the left with tufts attached to the trailing edge 
to verify that the flow remained attached and the tip vortex was steady. Vortex 
unsteadiness was a concern for the small airfoil aspect ratios (small values of 
penetration height, h) used to demonstrate the absence of blade vortex interaction 
noise when the vortex passed outside the rotor disk. 

Blade Tracking and Coning - Rotor blade tracking was evaluated during 
preliminary tests using a video camera and video monitoring screen. Synchronized 
light pulses simultaneously superimposed all four blades on the monitor screen. 

The camera viewing angle and synchronization provided a cross section view of the 
squared off blade tip. The sharp trailing edge on the monitor image indicated 
visually that all rotor blades were rotating in the same plane and blade pitch angles 
were identical. The accuracy of this visual assessment was based on observing a 
trailing edge thickness (of the superimposed blades) of less than 1 percent of the 
blade chord. This corresponds to a 16 minute variation between blade pitch angles 
which is within the - 10 minute limit used in setting pitch angle. The monitor 
image of the blade chord was 10 cm during this evaluation. 

Blade coning angles were also determined with the video system. Maximum 
coning angles of 0.5° indicated that the blades essentially rotated in a plane 
perpendicular to the rotor drive shaft. Although blade coning effects were small, 
they were accounted for when determining the vortex intersection angle, 3^» from 
flow visualization photographs. The procedure involved referencing 3^ to the local 
blade section in each photograph. 


Test Program 

Rational for Test Configuration - Before describing the rational for selecting 
the present vertical ascent test configuration it is necessary to describe the 
tail rotor blade vortex interaction noise mechanism. Figure 13a shows a simplified 
main rotor tip vortex trajectory as it is convected through the tail rotor. Tip 
vortices which initially are nearly normal to the tail .rotor plane of rotation 
(defined by the coordinates x^, X 2 for the case of zero cant angle on the rotor) 
are turned or skewed by the induced flow through the tail rotor. The skewed vor- 
tex velocity field subsequently creates an unsteady velocity component normal to 
the blade planform which results in noise radiation. Without turning of the flow, 
the vortex centerline remains normal to the plane of rotation and the transverse 
velcoity components are small. Under this condition tail rotor blade vortex inter- 
action noise is weak. (Note that small transverse velocity components can exist 
in this case, due to the blade pitch angle or twist and the axial mean velocity 
defect (if a defect exists) in the vortex core.) 
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The experimental method reported here permits varying the angle between the 
vortex centerline and the normal to the blade planform. The approach employs the 
simple but tractable test configuration shown in Figures 2 and 3. Here, skewing 
of the vortex centerline is achieved by contraction of the flow through the rotor 
as illustrated in Figure 13b. Vortex skewing relative to the blade normal now 
occurs in the plane defined by the X 2 , X 3 axes. Differences between the directions 
of the convected mean flow in Figures 13a and 13b only change the blade relative 
Mach number, and, hence, loading noise. The present test configuration, therefore, 
permits direct assessment of the noise mechanism. 

It should be noted that the model rotor shown in Figure 13b could be oriented 
with the rotor shaft at 90° to the open jet flow. This geometry would then simulate 
the tail rotor forward flight operating condition shown in Figure 13a. However, 
the upstream airfoil vortex trajectory would then be in the plane of the rotor 
instead of at 90° to the rotor as Figure 13a. Thus, the experimental arrangement 
in Figure 13b represents the closest simulation of the actual tail rotor inter- 
action with main rotor shed vortices. 

Coupling Effects Between Aerodynamic Parameters - It was desired to devise an 
experiment in which only one parameter varied. This was important because the 
vortex and rotor blade aerodynamic parameters were coupled. The following discus- 
sion defines the coupling effects which required* consideration when formulating 
the test matrix. 

First, changing the semispan airfoil angle of attack, o^, or the airfoil pene- 
tration height, h, changes the vortex strength. A further change in vortex strength 
would also appear to occur due to the mean flow contraction induced by the rotor. 
(This point is discussed in the section titled Velocity Field Incident on Rotor.) 

Second, changing blade number, B, pitch angle, a^, or tip Mach number, M^, 
changes the rotor induced contraction ratio. Consequently, vortex intersection 
station, I, and intersection angle, Bj> are changed. 

Due to the complicated coupling between aerodynamic parameters, experimental 
assessment of the blade vortex interaction noise mechanism necessitated a specific 
test program sequence. Parameters which had a minimum coupling effect were assessed 
first. For example, blade number was tested with a minimum change in vortex 
strength or intersection angle, Bi. Conclusions obtained from the blade number 
sensitivity test were then used in the assessment of other parameters. 

To illustrate the approach, blade tip Mach number dependence was investigated 
by changing M^ and blade number. The latter change was required to maintain a 
constant contraction ratio, thereby ensuring fixed values for the intersection 
station, I, and intersection angle, 3^. When necessary, airfoil penetration height, 
h, was also changed to hold I (or B^) constant. The subsequent change in vortex 
strength due to varying h was documented with the independent vortex velocity 
field measurements. 
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Formulation of Test Matrix - The range of the parameters investigated 
experimentally can now be specified and the rationale will be outlined where 
necessary. Open jet speed, U Q , was selected to be 9 m/sec. A low velocity was 
chosen since it created a large rotor inflow contraction ratio. This ensured 
that the vortex centerline would be skewed relative to the normal to the rotor 
blade planform. 

The effect of vortex strength was evaluated by setting different angles of 
attack on the upstream airfoil (<*2 - 0°, 6°, 12°). Flow visualization using tufts 
verified that the airfoil operated unstalled at all values of o^. Although NACA 
0012 airfoils typically operated unstalled for a 2 > 12°, larger angles of attack 
were not tested. This decision was based on the tuft motion which indicated the 
existence of an unsteady tip vortex at the low 9 m/sec open jet speed. Such 
unsteadiness could result in undesirable wandering of the vortex at the rotor 
plane. This condition was most noticeable at small airfoil penetration heights, 
h, corresponding to small airfoil aspect ratios. 

Blade vortex intersection station, I, was selected to be 1.25 rotor blade 
chords inboard from the blade tip. At this station, three dimensional tip effects 
not accounted for in the acoustic theory were expected to be weak. Intersections 
further inboard (larger values of I/R-^) were undesirable since the associated 
decrease in local Mach number results in weak blade vortex interaction noise rela- 
tive to blade broadband noise. More important, at larger values of I/R^ the 
rotor blade intersects the upstream airfoil two dimensional wake which introduces 
an additional noise mechanism. Blade intersection station was held constant by 
adjustment of blade penetration height, h, as described below in the discussions 
of the remaining parameters. 

Blade number dependence, B, was evaluated using a two and four bladed rotor 
configuration. Due to coupling between aerodynamic parameters, changes in B also 
changed the open jet contraction ratio and, hence, the vortex strength, intersec- 
tion station, I, and intersection angle, $]_• Changes in I and 3^ were countered by 
a small adjustment of the airfoil penetration height, h. Changes in vortex 
strength were shown to be negligible. 

The influence of vortex intersection angle, B-p was determined by combined 
changes in h, B, and a-^. Changes in h (or the separation distance, H, in Figure 7b) 
caused the vortex to be convected on different path lines into the rotor. For 
a pathline originating inside the rotor radius (H positive) the vortex was 
convected through the rotor with minimum contraction resulting in Bi =0°. In con- 
trast, displacing the airfoil to negative values of H in Figure 7b and changing from 
B = 2 to B = 4 created a higher contraction ratio. In this way 3-^ was varied from 
0° to 14.6°. Blade number effects were accounted for using the previous independent 
assessment of this parameter. Changes in vortex strength were shown to be 
negligible. 
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Blade pitch angle and Mach number investigations also required additional 
changes in B and h to hold I and constant. Pitch angle varied from 5° to 14° 
blade geometric angle of attack. Mach number varied over a range of 
M-p = 0.35 to 0.55. The lower value was limited by merger of the blade vortex 
interaction tones and the rotor broadband noise. The upper value was limited by 
the generation of blade thickness noise at higher Mach numbers. Introduction 
of a second competing noise mechanism would complicate efforts to isolate the 
blade vortex interaction noise mechanism. 

Forward Flight Effects - For tests conducted at freestream Mach numbers below 
0.1, measurements outside an open jet can be used to infer the acoustic source 
characteristics. This was verified experimentally by Schlinker and Amiet (ref. 26) 
using experiments designed to assess the effects of shear layer refraction, reflec- 
tion, and scattering of sound. Based on their conclusions, microphone corrections 
were not needed for the present experimental study. 


Qualification of Acoustic Test Procedure 

Several preliminary acoustic measurements were conducted to establish the 
quality and limitations of the acoustic measurements. Facility and rotor rig 
background noise are documented in Figure 14. Here the rotor rig was operated 
without a hub and the upstream semispan airfoil was removed. Rotor speed was 
varied over the selected test range with tunnel speed set at 9 m/sec. Achieving 
the low background noise shown in the present figure required wrapping the rotor 
rig with commercially available 0.15 cm lead sheets covered with sound absorbing 
foam. This covering (not shown in Figs. 2 or 3) attenuated noise generated 
by internal rotating components such as the bearings on the rotor rig. The 
resulting background noise level was sufficiently below the rotor broadband noise 
levels as indicated by acoustic measurements presented in later sections. 

The accuracy of setting blade pitch angle and the close blade tracking are 
documented in Figure 15a. Discrete tones shown here represent blade loading noise 
since the semispan airfoil was removed. Accurate blade-to-blade pitch angle 
settings and the close blade tracking is verified by the absence of a discrete 
tone at the rotor shaft frequency. Additional verification is evident in the 
absence of higher shaft harmonics. Finally, increasing blade number from B = 2 
to B = 4 (Fig. 15b) provided a 6 dB increase in loading noise. This increase 
occurs only if all blades are operating at the same pitch angles. As in 
Figure 15a, discrete tones at shaft harmonics are absent for B = 4. It should 
be noted that the steady loading tones decay rapidly with harmonic number as 
expected from theory. 
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Figure 15a also confirms the absence of higher harmonic narrowband random 
noise due to ingestion of the open jet turbulent shear layer at M T = 0.39. 
Increasing the blade tip speed to M T = 0.55, as shown in Figure 16, again verifies 
the absence of turbulence ingestion noise and demonstrates that close blade 
tracking was retained at high tip speeds. 
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VELOCITY FIELD INCIDENT ON THE ROTOR 


Objective 

The objective of the velocity measurements was to document the steady and 
unsteady flow field incident on the rotor. Steady velocities investigated 
consisted of the azimuthal velocity (Vq) and the axial velocity defect associated 
with the vortex. In addition, the two-dimensional semispan airfoil wake axial 
velocity defect was surveyed. Unsteady velocity measurements included the turbu- 
lence intensity and length scales in the viscous core. 

Measurements of the azimuthal velocity component were used to assess the vor- 
tex velocity field model in the blade vortex interaction theory. In particular, 
these measurements were needed to verify the radial dependence of the axisymmetric 
flow assumptions in Equation 22. In addition, the velocity measurements determined 
the maximum azimuthal velocity, Vq, at the edge of the vortex viscous core and the 
viscous core radius, r D . These parameters were experimental inputs in the assess- 
ment of the acoustic theory. 

In essence, the above described measurements of the azimuthal velocity field 
documented if the vortex was fully-developed. An experimentally generated fully- 
developed vortex was critical since the present noise theory does not treat the 
intermediate stage between the initially flat trailing vortex sheet and the fully 
rolled up vortex. Fortunately, experimental evidence indicates that vortex rollup 
occurs within a few blade chords. Thus, the vortex development was not expected 
to be a limitation in the present experiment since the rotor was stationed 5 chords 
downstream of semispan airfoil trailing edge. 

Documentation of the axial velocity defects in the vortex and two dimensional 
wake were required since these velocity fields could potentially introduce additional 
unsteady upwash components sensed by the rotor. The two dimensional wake was in- 
cluded in this evaluation since the rotor blade vortex intersection occurred at 
1.25 chords Inboard of the blade tip. Under this condition, it was possible that 
the rotor blade interacted with the two dimensional wake resulting in a wake chopping 
mechanism similar to that occurring in turbo-fan engines equipped with inlet guide 
vanes. Since the present vortex interaction theory models only the blade response 
to Vq, it was essential to document the various axial velocity defects. 

Measurements of the vortex viscous core turbulence intensity and length scale 
were needed to assess the importance of incident turbulence noise. Potentially, 
blade-to-blade correlations existed due to multiple chopping of turbulent eddies 
in the viscous core. Under this condition, narrowband random peaks could occur at 
the same blade passing harmonics generated by chopping the vortex azimuthal velocity 
field. 
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Approach 


Steady and unsteady velocity measurements were conducted at the rotor plane 
5 chords downstream of the semispan airfoil. All surveys were performed in a 
plane perpendicular to the open jet centerline with the rotor removed. 

Several reasons existed for removing the rotor during the velocity surveys. 
First, the present single element slanted hot film probe measurement technique 
assumes nonexistent radial velocities in the flow field. This condition was not 
satisfied with the rotor installed since the flow was contracted to skew the vor- 
tex centerline relative to the blade chord normal. A crossed-probe sensor can be 
used to circumvent this limitation but such a sensor averages radial velocity 
information (Appendix A) in the vortex core. A second reason for removing the 
rotor was the danger involved with measuring the vortex velocity field at the rotor 
face. During operation of the facility any mismatch between the rotor rotation speed 
and the tunnel speed results' in coning of the rotor blades. Potentially, the 
blades could strike the velocity sensing probe. A final reason for removal of the 
rotor was to avoid contaminating the unsteady hot film velocity measurements with 
the particle velocity generated by the rotor acoustic field. 


Vortex Azimuthal and Axial Mean Velocity Field 

Radial velocity dependence of the vortex azimuthal velocity component, Vq ,. is 
shown in Figure 17. Here, Vq, is normalized by the open jet freestream velocity, 

U 0> while radial distance from the vortex center (defined schematically in Figure 9) 
is normalized by the semispan airfoil. Measurements shown here were obtained in 
the horizontal plane defined by the x, z axes in Figure 9. The vortex center is 
defined as the point at which Vq is zero. This provides a check on the ability 
to independently identify the vortex center using the approach defined in Appendix A. 

The velocity profiles in Figure 17 correspond to two airfoil penetration 
heights, h, or airfoil to rotor tip separation distances, H. A high resolution 
radial traverse is provided for each value of H to show the velocity profile inside 
the vortex viscous core. In each case, airfoil angle of attack is constant and 
corresponds to a 2 = 12°. 

The profiles in Figure 17 show a systematic change In vortex structure over the 
range of separation distances tested. Decreasing the penetration height, or the 
airfoil aspect ratio (AR) , decreased the maximum tangential velocity. This is 
expected since AR controls the strength of the vortex sheet which rolls up into the 
tip vortex. Viscous core diameter, D, normalized by airfoil chord C2» is shown 
to remain approximately constant with D/c 2 = 0.09. 
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Precise symmetry of the velocity field does not exist since the maximum 
positive and negative values of Vq/U q differ. In addition, decay of the potential 
velocity field outside the viscous core is slightly different on the pressure and 
suction sides of the airfoil. The degree of asymmetry can be evaluated by comparing 
the measured velocity field with the analytical velocity field given by Equation 22. 
Figures 17a show the calculated Vq velocity profile based on the average value of 
the maximum positive and negative velocity and the measured core diameter (D/C 2 = 
0.091). Although complete symmetry does not exist the predicted velocity gradient 
in the viscous core and the potential flow decay rate were considered to adequately 
model the vortex. 

Further assessment of the azimuthal velocity field is provided by the multiple 
traverses shown in Figure 18. Here traverses were initiated at different radial 
stations designated by the numbers 1 through 8 in the small insert showing the 
airfoil. In each case, the traverse started at a station in the potential velocity 
field and penetrated the vortex core. The series of curves has been shifted to 
allow plotting all velocity profiles on the same page. The vortex center is indi- 
cated on each trace. All traces have been plotted showing the absolute value of Vq . 

Figure 18 indicates that the vortex is fully developed at the rotor plane since 
the velocity profiles indicate a symmetric flow. The only exceptions to this 
symmetry are the fine structure of the profiles and the non-zero values of Vq at 
large values of r/c 2 for radial traverses 6, 7 and 8. The non-zero azimuthal 
velocities are possibly due to the downwash created by the semispan airfoil. This 
explains the need for flow visualization to establish the vortex intersection angle, 
Y> in Figure 7. 

Velocity profile surveys such as those in Figures 17 and 18 were used to docu- 
ment the vortex viscous core diameter and the maximum azimuthal velocity. These 
results are summarized in Figure 19 with data points spanning the range of the air- 
foil angle of attack and aspect ratios tested in the acoustic study. Lines have 
been drawn through the data points to permit interpolating the data. This avoided 
performing the time consuming velocity profile measurements for every test condition 
in the acoustic study. 

The results in Figure 19 indicate that the normalized vortex core diameter 
varies between D/C 2 = 0.07 and 0.09. Expressed in terms of the rotor blade chord, c, 
these values would be D/c = 0.16 and 0.20. Thus, the vortex core is small compared 
to the rotor blade chord ensuring an impulse response in the acoustic signature. 

This shifts the tones in the blade vortex interaction spectrum to high frequencies 
which helps to distinguish this noise from loading noise. In this way, the blade 
vortex interaction mechanism can be isolated. 
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Velocity results in Figure 19 indicate that Vg/U 0 is a linear function of 
airfoil separation heights H or aspect ratio, AR. The only exception to this is the 
measurement at = 6° and AR = 3.1. Here Vq/U 0 falls below the straight line 
extrapolation. A similar nonlinear dependence occurs for the vortex core diameter at 
this test condition. These discrepancies are presently not understood. Blade vor- 
tex interaction noise data was not acquired at this test condition. 

Axial velocity dependence on distance from the vortex center is shown in 
Figure 20 for a range of airfoil separation heights, H. The profiles indicate a 
velocity defect in the viscous core with the defect increasing as H/c 2 changes from 
-0.55 to +0.88. However, the separation heights used in the assessment of the vor- 
tex interaction noise mechanism only varied from H/C 2 = 0 to +0.55. Over this 
range the axial velocity defect is approximately 5 percent of the mean velocity 
level. Note that the value of H/c 2 is given in Table 1 for each test condition. 

(One test was conducted at H/c 2 = -0.55 but this did not generate impulsive noise 
since the vortex intersected the rotor plane outside the rotor radius.) 

Values of U/U Q less than unity outside the vortex core in Figure 20 suggest a 
reduced open jet speed. However, at sufficiently large radial distances, U/U Q 
approached the freestream value. The nonuniform open jet velocity over the distance 
shown in Figure 20 was not expected to contribute to unsteady blade noise. 


Axial Velocity Defect in Two Dimensional Wake 

The upstream semispan airfoil wake velocity defect is shown in Figure 21. Here 
the normalized parameter, 1-U/U 0 , is presented as a function of transverse distance 
on a line perpendicular to the airfoil planform (based on a 2 = 0°). The measure- 
ment was conducted at 0.88 chords inboard of the airfoil tip. The transverse 
measurement distance is represented as a normalized radial distance to permit 
direct comparison with the vortex velocity field. 

Based on Figure 21, the maximum velocity defect associated with the two dimen- 
sional wake corresponds to a 15 percent decrease of the open jet speed. This exceeds 
the 5 percent axial velocity defect existing in the vortex viscous core. The airfoil 
wake is, therefore, potentially the strongest source of unsteady axial velocity 
defects incident on the rotor blade. Note that the location of the airfoil wake 
velocity defect is shifted relative to the vortex center, possibly due to the 
airfoil downwash. 

For comparison, the radial variation of the vortex azimuthal velocity component 
is also shown in Figure 21 for the same test condition. Here the magnitude of 
Vg/U 0 exceeds the wake velocity defect suggesting that the rotor blade unsteady 
upwash will be dominated by the vortex azimuthal velocity field. However, based on 
a following subsection, titled Unsteady Upwash Due to Blade Vortex Interaction, the 
fluctuating azimuthal velocity component sensed by the rotor blade is approximately 
VgsinBi w ^ ere is the blade vortex intersection angle. As an example, if 3^ = 14.6° 
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(which corresponds to Test Condition B to be discussed later) the normalized unsteady 
upwash velocity is (V 0 /U o ) sin 14.6° = 0.12. This value is now comparable to the two 
dimensional wake axial velocity defect and exemplifies the need to assess the noise 
generated by the rotor blades intersecting the wake velocity defect. 


Turbulence in Viscous Core 

Previous subsections concentrated on the mean velocity field associated with 
the vortex and the two dimensional wake. It now remains to document the turbulence 
in the vortex viscous core. Normalized turbulence intensities are presented in 
Figure 22 for the azimuthal (vq/U 0 ) and axial (u’/U Q ) velocity components. A com- 
parison of vq/U 0 and Vq/U 0 (Fig. 17) indicates that unsteady azimuthal velocities 
sensed by the rotor are dominated by Vq. A comparison of u’/U 0 and the mean 
axial velocity defect in vortex core and the airfoil wake indicates that unsteadi- 
ness due to these velocities is of the same order of magnitude. Thus, the axial 
turbulence component in the vortex viscous core in addition to the axial mean 
velocity defects must be considered in the assessment of potential noise sources 
contaminating the blade vortex interaction noise experiment. 

Autocorrelations of the axial component of turbulence were also obtained to 
determine the eddy length scale associated with the turbulence. Figure 23 shows 
a sample autocorrelation normalized by the zero delay time value. The Eulerian time 
scale represented by the area under Figure 23, is given by J ux = 7 msec. Invoking 
Taylor’s frozen flow hypothesis the axial integral length scale becomes A^ = U G J ux = 
6.3 cm. Based on approximately a 3 msec time between blade vortex intersections 
(corresponding to Test Condition B in Table 1) the eddy associated with the cal- 
culated length scale would be chopped twice. Thus, multiple intersections of the 
turbulence in the vortex core was considered a potential source of narrowband 
random noise at blade passing harmonics. 


Changes in Velocity Field Due to Contraction of Flow 

Vortex and axial wake velocity measurements were obtained with the rotor 
removed. The following discussion demonstrates how these measurements are applied 
in the presence of the' mean flow contraction created by the rotor. The discussion 
begins by considering the azimuthal velocity component, Vq, which forms the basis 
of the present acoustic theory. The affect of contraction on the axial velocity 
defect is described later. 
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The contraction skews the initially straight vortex in Figure 13b causing the 
vortex centerline to intersect the rotor at angles 3^ and y. The azimuthal velocity 
field, documented by hot film measurements (without contraction) in the plane normal 
to the vortex centerline, is similarly reoriented. The effect on the upwash velocity, 
w, at the rotor plane is accounted for theoretically by the intersection angles 
6 V and <t> v in Equations 36 and 38. Increases in the radial velocity field due to 
vortex stretching during the contraction are, however, not incorporated in the theory. 
These changes are known to occur based on Prandtl's model (ref. 27) for a contract- 
ing flow. His model predicts that the azimuthal velocity in an axial vortex increases 
as the square root of the velocity contraction ratio. Since this change could not 
be monitored experimentally it was necessary to analytically link the velocity 
fields occurring with the rotor removed and installed. 

Before developing the analytical correction to the measured velocity data, it 
is useful to examine in detail the changes in the unsteady upwash velocity due to 
changing the azimuthal velocity, Vq, or the vortex core radius, r Q . Based on 
Equation 36, changes in Vq or r Q only influence the parameters r Q V 0 and r Q ^. Thus, 
if the stretched vortex was also defined by the radial velocity function in Equa- 
tion 22, predicting the contraction effect would be limited to calculating the 
changes in r Q and V Q , 

It now remains to develop the equations defining the vortex after contraction. 
Assume initially that the radial velocity function in Equation 22 also models the 
stretched vortex. (The proposed similarity of the velocity function is a necessary 
condition for satisfying Kelvin’s theorem which is described below.) Then the 
parameters tq ^ and v 0,l Equation 22 are replaced by rg ^ and Vo,2* Here, 

the subscripts 1 and 2 are used to denote conditions before and after the contrac- 
tion. 


The product, r o 2 V 0 2 can to r o ,V 0 1 using Kelvin’s theorem 

(ref. 28). Based on this theorem, the circulation, T, around any contour line in 
the fluid remains constant in time as the contour is convected in an inviscid flow. 
Selecting a contour line at the edge of the viscous core and calculating the cir- 
culation before contraction gives T ^ = 277rQ^Vg Similarly, the circulation after 

contraction is = 2-rrrg 2 ^g 2* T ^ ie equality of T required by Kelvin's theorem 
results in rg 2 V o 2 = r 0 1 V 0 1* T Eus , the product r 0 V Q in the Equation for w remains 
invariant as the vortex is contracted. 

It should be noted that for an ideal inviscid potential vortex the product, 
rV, is invariant with radius and represents the vortex strength K^ = F/2tt. For 
this reason the product r 0 V Q is referred to as the vortex strength in the present 
study. Numerical values for the product are, however, expressed as a normalized 
vortex strength K 2 = ^ 702 ^. 
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The parameter r Q ^ in Equation 36 must also be linked with conditions prior to 
contraction. In this case, the incompressible continuity of flow condition can 
be applied. The radius of the vortex is then expected to decrease with the open 
jet contraction ratio given by r^/rg. Here, r-g and r^ represent conditions before 
and after contraction and the ratio is estimated using r^/r-g = (R^ - I)/(R-^ - H) . 
The vortex core radius r Q , therefore, varies as Tq 2 = r 0,l ^ r A^ r B^* 

The effect of contraction on the axial velocity defects in the vortex core 
and the two dimensional wake can also be predicted. Based on simplified continuity 
arguments the axial velocities increase in proportion to (rg/r^)^. 

In summary, the above discussions demonstrate that measurements of the vortex 
velocity field obtained with the rotor removed can be applied for the present con- 
tracted flow. Hot film measurements of the product r Q V 0 can be used directly as 
input to the non-contracted upwash velocity expression given by Equation 36. The 
effect of changing r Q due to contraction can be accounted for analytically before 
introducing the experimental value in Equation 36. 


Unsteady Upwash Due to Blade Vortex Interaction 

Before describing the experimental acoustic results of this study, it is useful 
to understand the unsteady upwash sensed by the rotor blade in the present vertical 
ascent test configuration. Figure 24 shows a plan view, front view and side view 
of a rotor blade during vortex interaction. (This figure is a detailed representa- 
tion of the blade vortex interaction for blade number 1 in Figure 7.) The vortex 
is represented by the azimuthal velocity field, V, shown in the plan view. The 
front view shows four vectors Vj_, V 2 , V 3 and which represent V at different 

azimuthal angles. These vectors lie in a plane perpendicular to the local vortex 
centerline and the radial distance from the vortex centerline is the same for 
each vector. Thus, the vectors are of equal magnitude. Note that the local blade 
angle of attack is zero to simplify this description, although the projected length 
of the vectors may not be equal. 

The side view of Figure 24 can be used to examine the unsteady upwash velocity 
components sensed by the rotor blade. Here the vortex, represented by the circula- 
tion vector, T, intersects the rotor blade at angle g-^. In the plan view the vor- 
tex arrives at y s 0 ° which corresponds to the present test conditions. 

For the zero angle of attack conditions shown here vectors V 2 and V 4 generate 
velocity components parallel to the blade chord. Vectors and V 3 , however, 
generate unsteady upwash components. For example, the side view in Figure 24 shows 
that the contribution due to Vj_ is - V-^sinf^. Based on this simplified figure, 
the unsteady upwash is a function of sing^. Each chordwise station experiences 
a simultaneous upwash and downwash due to vectors and Vg. Also, when the blade 
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is at angle of attack, vectors and contribute to the unsteady velocity field. 
These three dimensional effects illustrate the complicated unsteady flow associated 
with the blade vortex interaction mechanism and exemplify the need for the present 
theoretical formulation. 
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EXPERIMENTAL ASSESSMENT OF BLADE 
VORTEX INTERACTION NOISE 


Objective 

The objective of this phase of the study was to experimentally assess the 
noise generated by the blade-vortex interaction mechanism. Incident vortex charac- 
teristics and rotor operating conditions were varied to determine the sensitivity 
of the noise mechanism to isolated changes of each parameter. Parameters investi- 
gated included vortex strength, vortex intersection angle, blade number, blade 
pitch, tip Mach number, and directivity. Experimental assessment of each physical 
parameter involved conducting tests at two different operating conditions. Acous- 
tic spectra or acoustic pressure signatures were compared to determine the noise 
radiation sensitivity to the individual parameters. These experimental results 
also provided a direct assessment of the theoretical prediction procedure. 


Effect of Blade-Vortex Interaction 

Objective and Approach - The objective of this subsection is to describe the 
general features of blade vortex interaction noise. Acoustic spectra and pressure 
signatures are presented with the upstream vortex generator removed and installed. 
Measurements are also presented documenting the noise generated by the rotor blade 
intersecting the vortex axial velocity defect and the two dimensional wake defect. 
In addition, noise generation due to blade interaction with turbulence in the vor- 
tex core and two dimensional wake is assessed. Finally, the issues of spectrum 
analyzer filter bandwidth and the potential effect of small changes in rotor speed 
are addressed. 

Experimental Results - Figure 25 compares the acoustic spectra obtained with 
the upstream vortex generator removed and installed. For the installed case, the 
airfoil angle of attack was = 12°. The four bladed rotor tip Mach number was 
Mj, = 0.55 while the blade pitch angle was = 9.5°. (Additional parameters are 
listed in Table 1.) The microphone measurement angle in Figure 25 corresponds to 
a radiation angle, 0 O = 120°, in the downstream quadrant. This measurement angle 
was selected instead of the 0 Q = 90° station since the blade-vortex interaction 
dipcle radiation pattern has a minimum in the plane of rotation. 

Figure 25a shows that without the semispan airfoil the acoustic spectrum is 
dominated by the blade passing frequency (BPF) with higher harmonics decaying 
rapidly until they merge with the broadband background noise. Since the isolated 
rotor was operated at sufficiently low tip speeds to avoid generating high speed 
impulsive noise, the measured spectrum is considered to be dominated by loading 
noise . 
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With the semispan airfoil installed (Fig. 25b) the acoustic spectrum is 
significantly altered with the generation of numerous high frequency tones at 
multiples of blade passing frequency. Blade vortex interaction tones are 20 dB 
above the isolated rotor broadband noise level. Tones occur over a frequency 
range which permits distinguishing the interaction mechanism from the loading 
noise. 

A further understanding of this noise mechanism can be obtained by considering 
the acoustic pressure signatures. Figure 26 shows ensemble averaged pressure wave- 
forms with the semispan airfoil removed and installed. These time histories were 
obtained using the time marker signal from the rotor rig to initiate repetitive 
sampling of the radiated acoustic waveform. A 25 psec time increment between 
digital samples of the acoustic signature ensured adequate resolution of the 
waveform as the signal analyzer averaged 64 sequential time records. Polarity 
reversal of the acoustic signal due to the microphone signal conditioning elec- 
tronics was accounted for by reversing the standard sign convention when labeling 
the vertical axes in Figure 26. 

With the airfoil removed in Figure 26a, the sine wave like waveform represents 
the loading noise which dominated the spectrum in Figure 25a. Introducing the 
upstream airfoil (Fig. 26b) adds an impulsive signature to the measured waveform. 

Using the time marker from the rotor rig as a reference, the source of this impulsive 
noise was identified as the blade vortex interaction. This confirmed that higher 
harmonics in the spectrum of Figure 25b represent the vortex interaction mechanism. 

The above described differences between the loading and vortex interaction 
pressure signatures help to distinguish the two noise mechanisms. A further isola- 
tion is needed, however, if experimental results are to be compared with the present 
theory which models only the vortex interaction mechanism. Thus a procedure was 
developed to remove the loading noise contribution from the time waveform. 

The method is described in detail in Appendix B. 

Figures 26c and 26d show the resulting ensemble averaged blade vortex interaction 
noise signature for a single vortex intersection. (Note that smaller sampling incre- 
ments (5 ysec) were used for the high resolution time history in Figure 26d) . The 
acoustic response consists of a positive overpressure, followed by a negative over- 
pressure, and a second positive overpressure. The origin of this unsteady acoustic 
pressure signature was the angle of attack excursion sensed by each rotor blade 
as it intersected the vortex velocity field. As required by the acoustic noise 
model, the unsteady pressure is zero prior to and following the vortex intersection. 
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Importance of Axial Mean Velocity Defect - The above discussion does not 
address the potential noise generated by the rotor blades intersecting the axial 
velocity defects in the vortex viscous core and the two dimensional wake region. 

Flow measurements indicated these additional sources of unsteady upwash existed in 
the present experiment. The magnitude of their contribution to the measured 
acoustic spectrum therefore required assessment. This was achieved by repeating 
each acoustic measurement with the semispan airfoil at = 0°. Under this 
condition the axial velocity defects are approximately the same but the vortex 
azimuthal velocity component is zero. 

The spectrum in Figure 25 demonstrates that with the airfoil at = 0°, the 
higher harmonics generated by the blade vortex interaction are significantly weaker 
when compared to the a .2 — 12° spectrum. In general, the difference between = 0° 
and a 2 = 12° varied from 5 to 10 dB for all test conditions in the present study. 

This suggested that the measured acoustic response was dominated by the vortex 
azimuthal velocity component, Vg. It was, however, recognized that the airfoil 
wake defect could change as a .2 changed from 12° to 0° in the above described spectrum 
comparison. Thus, the two dimensional wake region could not be totally eliminated 
as a potential noise source. However, additional experimental results confirming 
the dominance of the vortex intersection over the two dimensional wake is presented 
in the subsection titled Vortex Intersection Angle Dependence. In this case, the 
upstream airfoil angle of attack (and consequently, the vortex strength) remained 
constant while the vortex intersection angle was changed at the rotor disk. 

It should be noted that the symbol for the ot 2 = 0° data in Figure 26 consists 
of a circle combined with a dashed line. The circle is used to identify the tone 
amplitudes while the dashed line defines the broadband noise level between tones. 

The absence of a dashed line indicates broadband levels remained essentially un- 
changed. For those regions of the acoustic spectrum where the decibel difference 
between the = 0° and = 12° spectra is small (for example, at 4 kHz in 
Figure 25b) experimental assessment of the theory was not attempted. 

Importance of Incident Turbulence Noise - Unsteady velocity field measurements 
inside the vortex core suggested that incident turbulence noise could, potentially, 
contribute to the measured acoustic spectra. Based on the length scale shown in 
Figure 23, eddies in the vortex could be chopped multiple times resulting in narrow- 
band random peaks at multiples of blade passing frequency. Similarly, turbulent 
eddies in the two dimensional wake could be chopped multiple times . Without an 
evaluation of this potential additional noise source the origin of the discrete 
tones in Figure 25b would be unknown. 

One method for identifying the presence of turbulence ingestion noise employs 
a high resolution narrowband frequency analysis. Using this approach randomness 
in the turbulent eddy structure is observed to cause broadening of the discrete 
tones generated by multiple intersections of the eddy. 
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Figure 27 shows a high resolution spectrum of selected tones appearing in 
Figure 25b. At large blade passing frequencies such as 7 x BPF or 13 x PBF, the 
apparent discrete tones in Figure 25b are shown to be narrowband random peaks. 

Although this observation suggests the existence of incident turbulence noise, 
several features must be carefully considered. First tone broadening due to 
multiple intersections with a turbulent eddy results in equal broadening at all 
blade passing harmonics as reported by Hanson (ref. 29). In contrast, the 
spectrum bandwidth in Figure 27 increases with frequency. Furthermore, Figure 
22 of reference 29 predicts approximately 30 eddy intersections are needed to 
generate a 10 Hz bandwidth at the 10 dB down point in Figure 27. This exceeds, 
by a factor of 10, the number of eddy intersections calculated from Figure 23. 

Finally, a 0.82 m axial length eddy would be required to generate the 30 inter- 
sections based on a freestream convection speed of 9 m/sec. Since the vortex 
core is onry 0 . 008 m in diameter , the eddy length-to-diameter ratio would be 
approximately 100 which is unrealistic. In summary, the physical constraints 
imposed by the eddy axial length scale invalidates incident turbulence as a 
significant noise source in the present study. Similar arguments and conclusions 
apply for the turbulence in the two dimensional wake. 

In the absence of turbulence ingestion noise it remains to establish the source 
of the narrowband random peaks in Figure 27. Noting the increase in spectrum band- 
width as tone frequency increased, it was conjectured that small variations in 
rotor speed broadened the tones generated by the vortex interaction mechanism. 

This possibility was evaluated by examining the spectrum of the rotor shaft time 
marker signal for similar broadening. Figure 28 shows the high resolution spectrum 
for the same center frequencies presented in the acoustic spectrum of Figure 27. 

For example, the 28th shaft harmonic corresponds to the 7th blade passing harmonic 
for the four bladed rotor. A comparison of Figures 27 and 28 shows that the spec- 
trum bandwidth, and even details of the spectrum shape, are the same. This establishes 
that small variations in the rotor speed (approximately ± 0.2 Hz based on Figure 28) 
created the narrowband random peaks observed in the acoustic spectrum. 

Importance of Filter Bandwidth in Comparison of Experimental Data With Theoreti- 
cal Predictions - Before the above described narrowband random measurements can be 
used for comparison with theory, the issue of filter bandwidth must be addressed. 

This is necessary since the present vortex interaction theory predicts only the 
acoustic energy generated at each blade passing harmonic. In essence, the theory 
predicts the total area under the distributed spectrum shape in Figure 27. Thus, 
comparison of measured and predicted data requires careful consideration. For- 
tunately, all data used for assessment of the present theory was obtained using a 
12.5 Hz bandwidth. This bandwidth encompasses the narrowband random peaks at all 
frequencies in Figure 25. Filter bandwidth, therefore, was not a limitation when 
comparing theory and experiment. 
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It is recognized that the tone broadening due to variations in the rotor 
rotation speed can be eliminated electronically. This requires using a commericial 
spectrum tracking unit in conjunction with the spectrum analyzer. Such a unit 
was not available during the acoustic data reduction phase of the present study. 

The above assessment of the filter bandwidth, however, establishes that the processed 
data can be used without introducing errors. 

The only exception to direct comparison of measured and predicted spectra 
arises when the narrowband tone (or even a discrete tone) spans two filter bands. 

This can occur irrespective of the use of a spectrum tracking unit. In this case, 
the measured power spectral density in adjacent bands must be added to give the 
total acoustic energy. Without this the tone amplitude can be 3 dB low if the tone 
is centered between filter bands. This may be the source of irregularities observed 
in the spectrum shape in the present study. 

E ffect of Rotor Speed Variations on Acoustic Signatures - In the absence of a 
spectrum tracking unit, the effect of a varying rotor rotation speed on the ensemble 
averaged acoustic signature must be assessed. For the 79 RPS test condition shown 
in Figure 26, the shaft frequency varied - 0.2 Hz. This variation is expected to 
occur over many revolutions since the rotor hub, shaft, and drive motor inertia 
damp the dynamic response of the system. 

Using a conservative assumption that the rotation speed change occurs over 
one revolution, the rotation period would be compressed or expanded by 32 psec. The 
impulsive signature in Figure 26c would then be shifted relative to the once per 
revolution time marker signal resulting in potential smearing of the ensemble averaged 
time waveform. Assuming the - 32 ysec change develops linearly over one revolution, 
this corresponds to an - 8 psec compression or expansion of each 1/4 period. Poten- 
tially, the first impulse in Figure 26c, which is enlarged in Figure 26d, will be 
shifted - 8 psec relative to the time marker signal. The effect is negligible com- 
pared to the 500 psec duration of the acoustic impulsive signature. 

Blade Tracking - Based on the similarity of the acoustic signatures in Figure 
26c each blade generates the same unsteady lift due to vortex intersection. Small 
differences between blades are controlled by the loading noise waveform as shown In 
Figure 26a. 


Blade Number Dependence 

Approach - Figure 29 shows the effect of blade number, B, on the blade vortex 
interaction noise spectrum. The upper curve corresponds to the four bladed rotor 
while the lower curve was obtained for the two bladed configuration. Since the 
mean flow contraction changed with B, an adjustment in airfoil separation distance, 
H, was necessary to maintain a constant intersection station, I. This adjustment 
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held I/R-^ constant while decreased from 14.6° to 11.9°. Vortex strength, r 0 V Q , 
changed only 3% based on the change in H while vortex radius, r c , remained essen- 
tially constant. All other parameters were identical as indicated in Table 1. 

Note that the above described vortex properties are based on velocity field 
measurements with the rotor removed. Earlier discussions of the vortex velocity 
field, however, established that the rotor induced mean flow contraction did not 
change the vortex strength, r 0 V 0 , at the rotor plane. Also, it was shown that 
the contraction effect on r Q can be calculated using the mean flow contraction 
ratio. (Resulting changes in acoustic radiation due to r Q changes were, however, 
shown to be negligible in the present study.) These are important conclusions 
which are applied throughout the following experimental assessment of the acoustic 
theory. All assessments used the calculated vortex properties at the rotor plane 
as experimental inputs to the theory. 

Results - The major effect observed in the comparisons of Figures 29a and 29b 
is a uniform decrease in noise level of approximately 7 dB as blade number changes 
from B = 4 to B = 2. To further quantify this result, differences between measured 
acoustic sound pressure levels are presented in Figure 29c at multiples of rotor 
shaft frequency. Here an average change of 7.5 dB is observed. 

Present theory predicts a 6 dB change at all frequencies as blade number 
changes from B = 4 to B = 2. The larger measured change in Figure 29c was traced 
to the 1.7° decrease of B]_ which also occurred as B changed. Including this 
decrease in 3-p the predicted change is given by the solid line in Figure 29c. 

The good agreement between theory and experiment is considered verification of 
the analysis considering the noise measurement accuracies and the variation in 
Bj_ . Figure 29c also verifies the well established 6 dB change associated with 
loading noise as blade number doubles. Experimental data at shaft harmonics 
given by mB = 4 and 8 represent the loading noise source. 

It is important to note that the observer input angle, 0 Q , in the theoretical 
prediction for Figure 29c was increased by the blade pitch angle, a^. The rational 
for this correction was discussed earlier in the section titled Description of the 
Experiment (Subsection - Acoustic Measurements) . This correction procedure is 
applied for all predictions presented in the remaining subsections. 


Vortex Strength Dependence 

Approach - Figure 30 shows the effect of vortex strength changes on blade vor- 
tex interaction noise. Here the upstream airfoil angle of attack was increased 
from (*2 = 6° to 12°. Measured vortex strength, r Q V 0 , was increased by a factor 
of two while vortex radius remained approximately constant. With the exception of 
the blade wake changes all other parameters were held constant. The wake changes, 
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however, had a minimal effect on the noise generation as confirmed by the vortex 
intersection angle measurements described in the following subsection. This 
is because the vortex azimuthal velocity component dominated the acoustic response 
during the experiment. 

Results - A comparison of Figures 30a and 30b indicates vortex interaction tone 
amplitudes increase as vortex strength increases. Sensitivity of the noise to 
changes in vortex strength is defined quantitatively in Figure 30c where an average 
increase of 5.5. dB is observed at frequencies above the loading noise (represented 
by mB = 4 and 8) . This increase agrees closely with the theoretically calculated 
change which is also presented in Figure 30c. The prediction is based on the 
measured values of vortex strength and vortex radius and is represented as a solid 
line in the figure. 

The measured change in Figure 30c confirms several parametric features of the 
present analysis. First the unsteady blade upwash given by Equation 36 is a 
linear function of vortex strength r 0 V 0 . Doubling vortex strength, therefore, 
increases the acoustic intensity by 6 dB when Equation 36 is combined with 
Equation 8. Second, changing r 0 V Q in Equation 36 changes the spectrum uniformly 
since wavenumber and vortex strength appear as isolated parameters. Finally the 
linear dependence on r Q V 0 and the uniform effect on the acoustic spectrum applies 
irrespective of the blade vortex intersection angles 0y or <£y. (This observation 
is based on examining Equations (36) and (40)). The above conclusions are important 
when considering changing the main rotor tip vortex properties or the vortex 
trajectory with the objective of reducing blade vortex interaction noise. These 
conclusions provide useful scaling procedures for guiding future designs. 


Vortex Intersection Angle Dependence 

Approach - Figures 31a and 31b compare measured acoustic spectra for different 
vortex intersection angles, 6^ = 14.6° and 0°. The change in 8-^ was achieved by 
decreasing blade number, B, and blade pitch angle, a^. These parametric changes 
eliminated the mean flow contraction, thereby, forcing 8^ to zero. A small change 
in separation height, H, was also required to hold the blade vortex intersection sta- 
tion constant. The upstream airfoil angle of attack, however, remained fixed so 
that the two-dimensional wake defect and vortex axial velocity defect were essen- 
tially the same in Figures 31a and 31b. 

The multiple parametric changes described above must be accounted for and are 
described here to document quantitatively the effect on Figure 31. First, the 
earlier assessment of blade number dependence demonstrated a 6 dB reduction in 
sound pressure level as blade number changed from B = 4 to 2. Second, changes 
in H caused a 9% increase of vortex strength in Figure 31b. Based on Equations 
36 and 8, this results in a 0.7 dB increase of sound pressure level irrespective 
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of changes in blade pitch- Finally, the 2.5° decrease in blade pitch angle between 
Figures 31a and 31b resulted in approximately a 1 dB reduction in sound pressure 
level. This is based on the independent assessment of pitch angle effects which 
showed a weak dependence on . In summary, blade vortex interaction tone amplitudes 
in Figure 31b should be increased by 6.3 dB to achieve the same test condition as in 
Figure 31a. This result will be used shortly to interpret the significance of 
Figure 31. 

A similar discussion of test conditions and changes in parameters is required 
before Figure 32 can be interpreted. Here, vortex intersection angle changed 
from = 14.6° to 3^_ = 4.4°. Similar to Figure 31, B, H, and were also changed. 
The acoustic responses to these changes essentially cancelled each other. 

The above described experimental changes in vortex intersection angle are 
based on smoke flow visualization photographs with the angle 3^ referenced to the 
rotor blade. In contrast, 0y is referenced to the normal to the blade planform. 
Angles 3-^ and 0y are equivalent only if a-^ = 0°. Since the present test program 
was conducted with finite blade pitch angles an analytical transformation was 
employed to convert the measured values for 3-^ and y to blade referenced angles 
0y and 4'^* 

Results - The 16 dB decrease due to changing 6^ between Figures 31a and 31b 
(after incorporating a 6.3 dB correction to account for other parameters changing) 
verifies a fundamental result of the present acoustic theory. The analysis pre- 
dicts that a vortex intersection normal to the blade planform does not generate 
noise. This is dramatically confirmed by the large reduction in sound pressure 
level in Figure 31b. The low sound pressure levels in Figure 31b also confirm 
that axial velocity defects did not significantly contribute to the tone noise 
generation in the present study. This is an important conclusion since it confirms 
that the vortex azimuthal velocity component dominated the acoustic response during 
the experiment. 

A quantitative assessment of the sensitivity to vortex intersection angle is 
obtained by comparing the tone amplitudes in Figures 32a and 32b. Note that a 
finite intersection angle was employed in Figure 32b to ensure that the vortex 
azimuthal velocity dominated over the axial velocity defects during the blade vor- 
tex interaction. Differences between tone amplitudes in these figures are pre- 
sented as a function of shaft harmonic number in Figure 32c. The difference has 
been corrected for blade number but changes due to varying H or a-^ are not included. 
Theoretical predictions (corrected only for blade number effect) are also shown in 
Figure 32. The predictions are based on the experimental test conditions and in- 
clude the measured vortex strength and a-^. The difference between the predictions 
for Figures 32a and 32b can, therefore, be compared directly with the experimentally 
determined differences in tone amplitudes shown in Figure 32c. The comparison of 
measured and predicted sensitivity due to changing 3^ shows reasonable agreement with 
theory. The observed change is considered verification of theory considering the 
accuracy of measuring 3]_* 
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It is interesting to note that the dB difference in Figure 32c approaches 
zero for shaft harmonic number mB = 4 and 8. This is because these harmonics 
are dominated by loading noise as indicated from comparing the - 12° and 0° 
in Figures 32a and 32b. 

Having confirmed the vortex intersection angle sensitivity, a simplified 
dependence on the angle 0y can be extracted from the general theory. Equation 
(36) indicates that when cos 0^ is approximately unity the unsteady upwash is 
controlled by the tan0y term. This dependence requires only that <f>y remain 
constant and applies irrespective ' of the absolute value of cj>y (a conclusion 
determined from combining Equations 40 and 36) . The dependence of tan0y is also 
uniform over all wavenumbers, or in other words, over the complete acoustic spectrum. 

To demonstrate the application of this scaling law consider estimating the 
change in Figure 32c. Initially, g^ = 14.6° and y = 0° corresponds to blade 
referenced intersection angles given by 0y = 17.4° and <J>y = 37.6°. As the inter- 
section angle changes to give = 4.4° and 0°, the blade referenced angles become 
0y = 8.3° and <J>y = 32.3°. The estimated change in tone sound pressure level 
associated with these different vortex intersection conditions is 6.6 dB . This 
value is comparable to the analytical prediction in Figure 32c. Note that the 
variation in <J> y could not be accounted for in the simple tan0y scaling law creating 
a difference between the predicted and estimated effect on the noise radiation. 


Blade Pitch Angle Dependence 

A pproach - Blade pitch angle dependence was determined by operating the rotor 
blades at two pitch angles given by = 9.5° and 14°. As in previous test cases, 
blade number and airfoil separation distance required adjustment to hold intersec- 
tion angle, gp, and intersection station, I, constant at the rotor plane. Angle 
g-L was held constant within 1.8 degrees while all other parameters remained unchanged. 

Results - A comparison of acoustic spectra in Figures 33a and 33b documents the 
effect of changing blade pitch angle. Quantitative assessment of the change is 
given in Figure 33c as a function of shaft harmonic number. A general observation 
is that blade vortex interaction noise is not significantly influenced by the small 
changes in a^. Figure 33c also shows the theoretically predicted change in tone 
level as pitch angle changes. (The prediction includes the small change in g^ to 
permit direct comparison with the experiment.) Based on the comparison with experi- 
ment, it is concluded that the theory adequately predicts the effect of pitch angle 
change. 

The analytically predicted effect of changing pitch angle can also be estimated 
from the previously identified tan0y parameter dependence. Changing by rotating 
the airfoil in Figure 4 through small angles about the y axis causes 0y to change 
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while remains approximately constant. Considering the intersection angles in 
Figure 33a, a-j_ = 9.5°, 3^ = 5°, and 7 = 0° correspond to 0y = 10.7° and <f> v = 27.9°. 
In Figure 33b, = 14°, 3^ = 6.8°, y = 0° corresponds to 0y = 15.5° and <}> v = 26.2°. 

Based on the simplified tan0 v dependence the increase in sound pressure level 
due to increasing ctj is estimated at 3.3 dB . This agrees with the experimental 
result and the predicted curve obtained from the computer coded version of the 
theory. Note also that changes in blade pitch angle influence the blade vortex 
interaction noise spectrum uniformly for small values of 0y. 


Mach Number Dependence 

A pproach - Sensitivity of blade vortex interaction noise to tip Mach number 
changes was assessed by varying M^ from 0.35 to 0.55. Blade number was also 
changed while the associated change in 3^ was limited to 1.2°. The remaining 
parameters were held constant. 

Results - The effect of increasing blade tip Mach number by 60% is demonstrated 
by comparing the measured spectra in Figures 34a and 34b. Differences between tone 
amplitudes in these figures are presented in Figure 34c. Note that data is pre- 
sented only at shaft harmonics for which a 5 dB difference existed between the 
a .2 = 12° and 0° specta. Recall that the olj = 0° spectra establish the contribution 
of the axial velocity defect to the radiated noise. 

The solid line in Figure 34c represents (including the small change in 6-^) the 
theoretically calculated change in sound pressure level due to increasing Mj. A 
4 dB difference between theory and experiment is observed at low frequencies but 
agreement improves at high frequencies. A decrease in the predicted curve with 
increasing frequency is attributed to the frequency dependence of the airfoil 
response function. 


Directivity Pattern 

Approach - Figure 35 shows measured vortex interaction acoustic spectra at 
three radiation angles. Spectra, such as these, were used to determine the source 
directivity pattern. Measurements were limited to a range of 0 Q = 60° to 140° 
due to the upstream and downstream anechoic chamber walls. In addition, the rotor 
rig support stand interferred with noise radiated to angles larger than 0 Q = 140°. 

All parameters were held constant in this evaluation. 

Directivity Pattern Based on Acoustic Spectra — For frequencies associated with 
blade vortex interaction noise (f > 3 x BPF) , tone noise spectra in Figure 35 exhibit 
a directivity pattern which is nearly symmetric about the rotor plane. A further 
evaluation of the vortex interaction noise directivity pattern is provided by 
Figure 36. Here measured tone levels are plotted as a function of 0 Q for increasing 
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blade passing frequencies. Loading noise dominated BPF, 2 x BPF, and 3 x BPF. 
Consequently, the directivity associated with the vortex interaction mechanism 
could not be investigated at these frequencies. However, at higher blade passing 
frequencies experimental evaluation was possible. In this case, the measured radia- 
tion patterns were found to be approximately symmetric about the plane of the 
rotor blade which corresponded to e Q = 80.5°. 

Predicted directivity patterns are also shown in Figure 36 for selected 
blade passing frequencies dominated by blade vortex interaction noise (5 x BPF, 

8 x BPF, 11 x BPF) . Good agreement between theory and experiment exists for the 
directivity pattern shape. In addition, the measured and predicted minima occur 
near the plane of the blade. Since the acoustic source dipole axis is aligned 
normal to the blade planform in the present analysis, the source radiation is zero 
in the plane of the blade as verified by the experiment. 

Although generally good agreement between experiment and theory was observed 
for directivity pattern shape large differences occurred in the absolute tone 
amplitude. These differences were due to a general decrease in predicted tone 
amplitude with increasing frequency for all radiation angles, 0 . This observa- 
tion is based on a comparison of measured and predicted acoustic spectra which is 
discussed later in the subsection titled Assessment of Acoustic Spectrum and Pres- 
sure Signature. 

It should be noted that the value of 0 q used for the input to the analytical 
prediction was increased by the value of a-^. This is due to the difference in the 
reference axes selected for the theoretical and experimental definition of 0 Q (see 
subsection titled Acoustic Measurements). According to Figure 4, 0 Q in the 
analysis is measured from the normal to the blade planform. In contrast, the 
experimental program referenced 0 O to the rotor centerline. At the blade vortex 
intersection station, these angles differed by the blade pitch angle setting. Com- 
parison between theory and experiment, therefore, required increasing 0 Q by 
ct-^ = 9.5° in the predictions. This approach is used consistently throughout the 
report . 

Directional Radiation Characteristics of Blade Vortex Interaction Impulsive 
Signature - The previous discussion of directivity pattern focused on tone amplitude 
changes occurring over the range of 0 Q = 60° to 140°. Additional directivity 
features of blade vortex interaction noise can be identified from acoustic pressure 
signatures. Figure 37 shows ensemble averaged waveforms measured at selected 
radiation angles 0 O — 60°, 90°, 120° and 140°. The data corresponds to the same 
test condition described in Figure 36. 

Waveforms in Figure 37 contain the combined loading and blade vortex interaction 
signatures. Dominance of loading noise is evident in the sine wave character of the 
time history at all radiation angles. In contrast the impulsive signature peaks 
at angles outside the plane of rotation but diminishes to zero at 0 O = 90° which is 
close to the plane of the rotor blade. These time history traces provide additional 
explanations for the directivity features determined from acoustic spectra. 
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Blade vortex interaction acoustic signatures are further isolated in Figure 38. 
Here measured waveforms are shown for the same microphone angles shown in Figure 37. 
Loading noise has been removed using the procedure described in Appendix B. Several 
waveform features are noteworthy. First, waveform shape is observed to be anti- 
symmetric about the plane of rotation. Second, peak amplitude decreases near the 
plane of rotation. Finally, absolute waveform amplitudes are equal on opposite 
sides of the rotor. These important characteristics confirm the dipole acoustic 
source model developed in the present study. 


Assessment of Acoustic Spectrum and Pressure Signature 

Obj ective - Previous discussions focused on the experimental assessment of the 
sensitivity of blade vortex interaction noise to isolated changes of each parameter 
in the acoustic theory. Conclusions obtained from this assessment confirmed that the 
theory models the parametric dependence of each individual parameter. It now 
remains to determine if the analysis predicts the measured acoustic spectrum shape 
and absolute sound pressure levels. In addition, the analysis will be tested to 
determine if the measured acoustic pressure signatures can be predicted. 

Acoustic Spectrum - Figure 39 shows measured and predicted spectra for Test 
Condition B. Before comparing these spectra, the presence of steady loading noise 
must be discussed since this mechanism dominated frequencies given by BPF, 2 x BPF, 
and 3 x BPF. This additional noise source was removed from the measured spectrum 
by first using the waveform subtraction technique described in Appendix B.l. A 
Fast Fourier Transform (FFT) of the resulting isolated blade vortex interaction 
time history (consisting of 13 consecutive vortex interactions in the continuous 
time waveform) provided the acoustic spectrum shown in Figure 39. This spectrum 
coincides with the spectrum in Figure 25b except for the BPF tone and the first 
two harmonics . 

The general conclusion from the spectrum comparison in Figure 39 is that 
tone amplitudes are predicted to within 5 dB at low frequencies. At high fre- 
quencies, the present analysis significantly underpredicts the data. 

Acoustic Signatures - Figure 40a shows the predicted acoustic signature for 
the experimental test case presented in Figure 38c. Both waveforms are plotted on 
the same scales and correspond to the same measurement angles. A comparison of 
the measured and predicted waveforms indicates that the duration of the acoustic 
signature is similar as is the sign of the maximum peak pressure. In addition, 
the experimentally observed antisymmetric waveform occurring on opposite sides of 
the rotor are predicted by the present analysis. Calculated waveforms were found 
to be identical (not shown) at 0 q = 50° and 120°, except for the sign reversal in 
the acoustic pressure. (Note that these angles are symmetric around the plane of 
the rotor blade in the present experiment.) The qualitative agreement between 
theory and experiment again conforms the dipole acoustic source model. 
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While the general features of the blade vortex interaction signature are 
predicted, details cannot presently be calculated. For example, the peak ampli- 
tude, triangular shape, and combined positive and negative peaks of Figure 38c 
are not predicted. 

Assessment of Dif f erence Between Theory and Experiment - The rotor blade inter- 
action with the two dimensional wake of the upstream airfoil was initially con- 
sidered to be a possible source of the difference between theory and experiment. 

This mechanism was, however, evaluated experimentally during the blade vortex 
intersection angle tests. The results indicated that when (3^=0, the axial 
velocity defects did not significantly contribute to the tone noise generation. 

This conclusion confirmed that the vortex azimuthal velocity component dominated 
the acoustic response during the present study. 

A second possible source of the difference between theory and experiment is 
the infinite span airfoil assumption used in the analysis. With this assumption, 
physical constraints imposed by the finite span rotor blade are not satisfied by 
the theory. To understand the consequences of the infinite span airfoil model, 
recall that the far field sound is related to the time derivative of the airfoil 
loading. Calculating the area under the far field pressure waveform is, therefore, 
equivalent to a time integral of the pressure signature in Figure 40a. The 
result of the integral is the difference between the lift when the airfoil is at an 
infinite distance upstream and downstream of the airfoil 


Ldt = L - L (41) 

CO —00 

For the present infinite span airfoil model and the 1/r decay of the azimuthal 
velocity vector, Equation 42 gives a finite or non-zero difference between L^, 
and L^. This occurs because the calculated pressure pulse in Figure 40 is always 
less than zero. Pressure peaks of opposite sign, as in the experiment, would 
cause the area in Equation 41 to approach zero. 


AREA = 


00 

/ 


The non-zero value obtained from Equation 42 represents a theory limitation 
since a finite span airfoil situated at an infinite distance from the vortex has 
a zero value for the quantity - L^. To approximate the finite span blade 
tested in the present experiment. Equation 22 was multiplied by an exponential 
decay function defined as exp [- (R 0 /r) 2 ] . Here R Q represents the radius at which the 
exponential amplitude is e~l. The rapid decay of the vortex velocity field due to 
the additional exponential decay function ensures that the integral in Equation 42 
has a zero value. Introducing the exponential decay function represents more than 
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a mathematical termination of the influence of the vortex velocity field. It is also 
an approximate simulation of a finite span airfoil. The acoustic signature is the 
same whether the vortex velocity field is terminated or the airfoil is terminated 
as in the case of a finite span. This equivalence, however, assumes that the blade 
end effects can be ignored. 

Using the exponential decay function to model the finite span blade results in 
predicted acoustic waveforms which display some of the features of the measured 
waveforms. This is illustrated in Figure 40b where predicted pressure signals 
are presented for a range of R Q /c values where c corresponds to the blade chord. 

Here R Q /c = 1 approximates the distance, I, in the present experiment, from the 
vortex interaction station to the blade tip. Terminating the analytical infinite 
span airfoil model at this distance generated an acoustic pressure signal which 
exhibits combined positive and negative peaks similar to the experiment. The 
r / c = 1 curve displays an initial positive pressure pulse, a negative peak, 
and a small positive pressure waveform which returns to zero as time increases. 
Increasing R c /c causes the acoustic pressure to return to the one-sided waveform 
.predicted in Figure 40a. 

The above discussion illustrated how the waveform shape depends on the infinite 
span airfoil assumption. This assumption can also result in incorrect predictions 
of the amplitude of the acoustic signature as illustrated by the following example. 
Consider a plane normal to the span line of the airfoil. Also let the normal to 
the airfoil planform and the vorticity vector be in this normal plane (<£ v = 0° 
or 180° and 8y arbitrary) . Then the unsteady transverse flow sensed by the airfoil 
will be antisymmetric about this normal plane; i.e., if the flow is upward on one 
side of this normal plane it is downward on the opposite side with an equal magni- 
tude. An observer stationed anywhere in the normal plane will measure zero sound 
since a positive pressure perturbation produced on one side of the normal plane is 
cancelled by the opposing negative pressure. 

If the infinite span airfoil in the above description is now truncated on one 
side of the normal plane (as in the case of the finite span rotor blade used in the 
present experiment) the opposing unsteady pressure responses no longer occur and 
the above described acoustic cancellation disappears. This finite span effect 
is not accounted for by the previously described experimental decay factor since 
the factor retains the symmetry of pressure response on the airfoil surface. Thus, 
it is possible to also achieve an increase in predicted sound level when using 
a truncated (finite span) airfoil. 

The asymmetry associated with the finite span provides a potential explanation 
for the discrepancy between predicted and measured acoustic spectra (Fig. 39) 
in addition to acoustic signatures (Fig. 40). This follows since the present 
microphone measurement plane almost coincided with the above described normal 
plane. Also, for the small values of (3^ tested in the present experiment, the 
vorticity vector was almost in the normal plane. 
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CONCLUSIONS 


1. Ingestion of a vortex by a rotor generates harmonic noise which extends to 
high frequencies. Acoustic pressure signatures associated with this noise mechanism 
display impulsive waveforms superimposed on the blade loading acoustic pressure 
signals . 

2. Linear dependence of measured blade vortex interaction noise on blade number 
and vortex strength was closely predicted by the present first principles acoustic 
theory. The absence of noise for vortex trajectories normal to the blade was also 
predicted. For oblique intersections with a blade, calculated changes in acoustic 
radiation due to intersection angle modifications were within 2 decibels of the 
measured data. A simplified scaling law was determined for estimating these changes 
at small intersection angles. 

3. The predicted low sensitivity to blade pitch angle modifications for a typical 
range of rotor operating conditions was confirmed experimentally. The strong noise 
dependence on local Mach number was calculated although the sensitivity to Mach num- 
ber change was overpredicted at low frequencies. 

4. Measured blade vortex interaction directivity patterns confirmed the theoretical 
dipole acoustic source model at all frequencies. The minimum sound pressure level 
observed experimentally in the plane of the blade was predicted analytically. 

5. Absolute levels for acoustic spectra and directivity patterns were predicted, 
without the use of empirical or adjustable constants, to within 5 decibels at low 
frequencies. At high frequencies, sound pressure levels are currently underpredicted. 

6. General features of the measured blade vortex interaction signature were pre- 
dicted. Antisymmetry of the measured waveform on opposite sides of the rotor plane 
confirmed the theoretical dipole directivity model. Duration of the impulsive 
acoustic signature in addition to the sign of the maximum pressure peak was predicted. 
Detailed features of the measured acoustic signature, such as peak amplitude, tri- 
angular shape, and combined positive and negative peaks, are not presently calculated. 
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APPENDIX A 


Vortex Velocity Field Measurement Technique 

Decomposition of Total Velocity Vector - The azimuthal (V 0 ) and axial (U) 
velocity components of the vortex velocity field are shown in Figure 9. Decomposi- 
tion of the measured voltage from the single slanted hot film into Vg and U 
employs the following procedure. 

First, it is necessary to recognize that the probe sensing element can be ro- 
tated until it lies in the plane formed by the components V 0 and U. (Note that 
the vortex radial velocity is assumed to be zero in this velocity decomposition 
approach.) This unique orientation is illustrated in Figure A-l. Here the V 0 and 
U components lie in the plane formed by the hot film sensor and the support prongs. 

In Figure 9 this corresponds to alignment of the probe sensing plane with a plane 
parallel to the z, y axes. In this case, the probe sensing plane is normal to a 
radial line penetrating the vortex center. 

To decompose the total velocity vector, Uj, in Figure A-l into the corresponding 
values of V 0 requires two separate measurements. These correspond to two orienta- 
tions which are 180° apart in Figure A-l. The velocity components are then given 
by the equations: 

V 0 = Urj, cos (3p-a-^) = U,j, (cos 3 ^ cos a^ + sin 3 ^ sin a^) (A-l) 

U = U T sin ( 6 F -a 1 ) = U T (sin cos a-j^ - cos 3 F sin a-^) (A-2) 

Representing the two separate normal velocity components at each wire orienta- 
tion as and gives 



U N1 = U T COS a l 

(A-3) 


U N2 " U T S±n a l 

(A-4) 

Since the sensor 
linearly related 

mean voltage output E-^ and E£ at each of 
to the normal velocity component 

the two orientations is 


u N i - KE i 

(A- 5) 
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where K represents the velocity versus voltage calibration coefficient. 
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Using Equations A-5 and A-6 in the expressions for Vg and U gives the sum and 
difference expressions for velocities measured using a standard crosswire 

V Q = 0.707 K(E 1 ~E 2 ) 

U = 0.707 K(E 1 +E 2 ) 

Thus, a measurement at two orientations in a plane normal to a radial line pene- 
trating the vortex center gives Vq and U. 

Instead of rotating the sensor at each radial velocity station to obtain Vg 
and U the following procedure was used. A continuous traverse was conducted on 
a radial line through vortex center for each of the two orientations. An analog 
output from the traverse system was used to drive one axis of a two directional 
plotter. Two continuous traces of the voltages E^ and E£ were then obtained as the 
sensor penetrated the vortex. Graphically summing and differencing the plotter 
traces decomposed the measured voltages into E^-E 2 and E 1 +E 2 at each radial sta- 
tion. A simple application of the calibration factor K then permitted changing 
the plotter axes notation to Vq and U as shown in Figure 17. 

Identification of the Vortex Center - The above described radial velocity 
traverse assumes a knowledge of the line penetrating the vortex center. The follow- 
ing simple procedure was developed to identify the vortex center. 

As described in the previous subsection, there exists an alignment of the 
sensor measuring plane for which Vg and U are easily determined. This plane is 
normal to the radial line through the vortex center. The plane can be identified 
by rotating the sensor 360° about the probe shaft. When the sensor voltage output 
reaches a maximum, all velocity components lie in the plane of the sensor. This 
orientation is designated in Figure A-l as the maximum voltage output orientation. 
Other orientations can be shown to give smaller voltage values with a minimum 
occurring at 180° from the maximum orientation. 

The vortex center then lies on a line perpendicular to the unique measurement 
plane. Rotating the sensor at numerous stations on the perimeter of the vortex 
provided a series of radial lines whose intersection point defines the vortex center. 
Figure A-2 provides a schematic description of this experimental procedure. A 
digital angular position encoder with a 0.33° resolution was used to sense the 
probe orientation. 
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APPENDIX B 


Experimental Technique for Isolating Blade^Vortex 
Interaction Impulsive Signatures 


Measured acoustic pressure signatures contained both loading noise and the 
blade- vortex interaction impulsive signature. The following procedure was employed 
to isolate the impulsive waveform permitting direct comparison of measured and 
predicted blade-vortex interaction pressure signatures. 

The experimental technique begins with the ensemble averaged waveform for the 
combined loading and vortex interaction noise. A sample time history is shown in 
Figure B-la where the time marker signal from the rotor rig is used to initiate 
repetitive sampling of the acoustic waveform. This waveform is obtained with the 
upstream semispan airfoil at a finite angle of attack (012 = 12° in the present 
example) . 

A second waveform is obtained with the airfoil at = 0° (Fig. B-lb) . Under 
this condition, the vortex strength is forced to zero while the loading noise 
remains the same. At the same time, noise generated by intersecting the axial 
velocity defect in the vortex core and in the two dimensional wake are considered 
to remain the same. The difference between Figures B-la and B-lb represents the 
isolated blade vortex interaction signature which is shown in Figure B-ld. An 
enlarged version of this waveform is presented in Figure B-le. The above described 
ensemble averaged pressure signatures and the arithmetic subtraction of these wave- 
forms were obtained on a dual channel Fast Fourier Transform spectrum analyzer with 
time domain averaging capability. 

The rationale for using the = 0° waveform is based on the need to account 
for the axial velocity defect which is also present in the waveform of Figure B-la. 
The noise generated by this velocity defect creates the small pressure pulse 
observed in Figure B-lb. Using the loading noise pressure signature with the 
upstream airfoil removed (Fig. B-lc) would fail to account for this noise 
mechanism. 
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EXPERIMENTAL TEST CONDITIONS 
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WEAK 

STRONG MAIN MAIN ROTOR 

ROTOR BLADE BLADE VORTEX 


FLIGHT DIRECTION VORTEX INTERACTION INTERACTION 



a) ILLUSTRATION OF MAIN ROTOR AND TAIL ROTOR BLADE VORTEX INTERACTION 


VORTEX 



b) ILLUSTRATION OF MAIN ROTOR BLADE VORTEX INTERACTION WITH 
VORTEX PARALLEL TO LEADING EDGE 


Figure 1 — Helicopter Blade Vortex Interactions 
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Figure 3 — Upstream Model Airfoil and Rotor Blade with Boundary Layer Trip 
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Figure 4 — Airfoil-Vortex Geometry 




SIDE VIEW 


Figure 5 — UTRC Acoustic Research Tunnel 





Figure 6 — Acoustic Research Tunnel Control Room with Data Acquisition 
and Data Reduction Instrumentation 
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b) SIDE VIEW SHOWING OPEN JET TEST FACILITY 


Figure 7 — Detailed Schematic of Vortex Generator and Model Rotor 
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UPSTREAM AIRFOIL 
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Figure 8 — Microphone Arrangement 
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2) FIGURE ASSUMES DOWNWASH DUE TO 
AIROFIL NON-EXISTENT 


Figure 9 — Definition of Vortex Velocity Field 
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HOT WIRE 
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Figure 10 - Tip Vortex Generating Airfoil and Probe Traverse System in Acoustic Tunnel 





Figure 11 — Smoke Flow Visualization Showing Tip Vortex for 
Two Operating Conditions 
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Figure 12 — Flow Visualization of Blade Vortex Interaction for Lifting Rotor 
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b) VORTEX TRAJECTORY IN OPEN JET WIND TUNNEL 
WITH ROTOR 

Figure 13 — Equivalence of Tail Rotor Vortex Interaction Process for Full Scale 
Helicopter and Open Jet Acoustic Wind Tunnel Model Rotor Experiment 


73 



TEST CONDITION: 
NO AIRFOIL 
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Hub or Blades 




TEST CONDITION: 
NO AIRFOIL 
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Figure 14 — Concluded 




TEST CONDITION: 
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a) B = 2 BLADES 

Figure 15 — Acoustic Spectrum Demonstrating Blade Tracking Accuracy and Absence of 
Freestream Turbulence 




TEST CONDITION: IDENTICAL 
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Figure 16 — Acoustic Spectrum Demonstrating Blade Tracking Accuracy and Absence of 
Freestream Turbulence 




TEST CONDITION — IDENTICAL 
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Figure 16 — Concluded 
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Figure 17 — Continued 
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NORMALIZED VELOCITY, V„/U 
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NORMALIZED RADIAL DISTANCE, r/c 2 
c) H/c 2 = 0.55 


Figure 20 — Axial Velocity Dependence on Radial Distance from Vortex Center 
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Figure 20 — Concluded 
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NORMALIZED VELOCITIES 



-°- 25 0 0.25 0.5 0.75 1.0 

NORMALIZED RADIAL DISTANCE, r/c 2 

Figure 21 — Comparison of Vortex and Two Dimensional Wake Velocity Field, H/c 2 = 0.55, 
AR =2.77, a ? = 12° 
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Figure 24 — Schematic of Vortex Velocity Components 
Responsible for Noise Generation 
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SPL IN dB AT 12.5 Hz BANDWIDTH 



Figure 25 — Effect of Vortex Generator on Acoustic Spectrum 
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a) VORTEX GENERATOR ABSENT 
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b) VORTEX GENERATOR INSTALLED AT c* 2 = 12° 


Figure 26 — Effect of Vortex Generator on Acoustic Signature 
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ACOUSTIC PRESSURE IN 




d) HIGH RESOLUTION TIME HISTORY OF BLADE VORTEX INTERACTION 

Figure 26 — Concluded 
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Figure 27 — High Resolution Spectrum of Narrowband Random Tones in 
Acoustic Spectrum, Test Condition B, 0 o = t2O°, 0 o = l8O°, 
Test Condition B 


95 



RELATIVE VOLTAGE LEVEL IN dB AT 0.2 Hz BANDWIDTH 




a) SHAFT FUNDAMENTAL b) 28X SHAFT FUNDAMENTAL 



Figure 28 — High Resolution Spectrum of Shaft Rotation 
Time Marker Signal, Test Condition B 
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Figure 29 — Concluded 
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Figure 30 — Concluded 
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mB, m = INTEGER 


c) DIFFERENCE BETWEEN SHAFT HARMONICS FOR 
INTERSECTION ANGLES /3 1 =14.6° AND 4.4° 

Figure 32 — Concluded 


103 


SPL IN dB AT 12.5 Hz BANDWIDTH 


70 


TEST CONDITION: H 
0 o = 6O°, <t> Q = 120° 



55 


45 


35 


0 1 2 3 4 5 

FREQUENCY, KHz 
b) a-j =14° 


Figure 33 — Blade Pitch Angle Dependence 
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Figure 33 — Concluded 


105 










DIFFERENCE, dB 


NOTATION 


Q EXPERIMENT 
THEORY 



mB, m = INTEGER 

c) DIFFERENCE BETWEEN SHAFT HARMONICS FOR 
MACH NUMBERS My = 0.35 AND 0.55 


Figure 34 — Concluded 
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Figure 35 — Comparison of Spectra at Different Measurement Angles 
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Figure 35 — Concluded 
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SPL IN dB AT 12.5 Hz BANDWIDTH 
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b) 2 x BPF 

Figure 36 — Tone Directivity Pattern 
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Figure 36 — Continued 
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Figure 37 — Concluded 
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Figure 38 — Radiation Pattern of Blade Vortex Interaction Acoustic Signature 
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Figure 39 — Comparison of Measured and Predicted 
Tone Sound Pressure Levels 
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Figure 40 — Predicted Blade Vortex Interaction Acoustic Signature 
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MINIMUM VOLTAGE ORIENTATION (ASSUMES ZERO RADIAL VELOCITY) 

NOTE — HOT FILM SENSING ELEMENT LIES IN A PLANE PARALLEL TO z. y PLANE IN FIGURE 9 

Figure A-1 — Velocity Components Sensed by Single Slanted Hot Film Probe 
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Figure A-2 — Schematic of Procedure Used to identify Vortex Center 
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b) ACOUSTIC SIGNATURE FOR STEADY LOADING AND BLADE 
VORTEX INTERACTION, a 2 = 0° 

Figure B-1 — Sequence ol Acoustic Signatures Showing Procedure for Isolating Blade Vortex 
Interacton Time History 
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Figure B-1— Concluded 
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